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ABSTRACT 

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev-Zeldovich (SZ) effect detections using the first 
15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) 
sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly 
through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we 
divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest 
all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) x 10^^ M©. Confirmation of 
cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of 
the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This 
gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed 
SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all 
the Planck clusters not previously detected in X-ray surveys. 

Key words, large-scale structure of Universe - Galaxies: clusters: general - Catalogs 
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1. Introduction 

This paper, one of a set associated with the 2013 release of data 
from the Planck^ mission (Planck Collaboration I 2013), de- 
scribes the construction and properties of the Planck catalogue 
of SZ sources (PSZ). 

Clusters of galaxies play several important roles in astro- 
physics and cosmology. As rare objects, their number density is 
especially sensitive to properties of the cosmological model such 
as the amplitude of primordial density perturbations (Peebles 
1980), and their development with redshift probes the growth of 
cosmic structure, hence perhaps helping to distinguish between 
dark energy and modified gravity explanations for cosmic accel- 
eration (e.g., see reviews by Borgani & Kravtsov 2009; Allen 
et al. 2011). The galaxies, hot gas and dark matter held in their 
gravitational potential wells provide a sample of the universal 
abundance of these components (e.g., Voit 2005), while the ther- 
mal state of the gas probes both the cluster formation mecha- 
nism and physical processes within the cluster such as cooling 
and energy-injection feedback (e.g., reviews by Fabian 2012; 
McNamara & Nulsen 2012). The study of the constituent galax- 
ies, including the brightest cluster galaxies normally found at 
their centres, allows sensitive tests of galaxy formation models. 

Because of these uses, there is considerable interest in devel- 
oping large galaxy cluster catalogues that can be used for pop- 
ulation and cosmological studies (e.g., Schuecker et al. 2003; 
Bohringer et al. 2004). Clusters are genuinely multi- wavelength 
objects that can be selected in several ways: optical/infrared 
(IR) imaging of the galaxy populations; X-ray imaging of 
bremsstrahlung radiation from the hot cluster gas; and through 
the Sunyaev-Zeldovich (SZ) effect (Sunyaev & Zeldovich 1972; 
Sunyaev & Zeldovich 1980) whereby scattering of cosmic mi- 
crowave background (CMB) photons from that hot gas distorts 
the spectral shape of the CMB along lines of sight through clus- 
ters and groups. 

Construction of cluster catalogues in the optical/IR and in 
the X-ray are relatively mature activities. The first large opti- 
cal cluster survey is now over 50 years old (Abell 1958; Abell 
et al. 1989), and current catalogues constructed from the Sloan 
Digital Sky Survey data contain over a hundred thousand clus- 
ters (e.g., Koester et al. 2007; Wen et al. 2012). In X-rays, 
large samples first became available via ROSAT satellite ob- 
servations (e.g., Vikhlinin et al. 1998; Bohringer et al. 2000; 
Gioia et al. 2003; Bohringer et al. 2004; Burenin et al. 2007; 
Ebeling et al. 2007), but also more recently for instance from 
dedicated or serendipitous survey with XMM-Newton (Pacaud 
et al. 2007; Fassbender et al. 2011; Takey et al. 2011; Mehrtens 
et al. 2012). Currently several thousand X-ray selected clus- 
ters are known (see for instance the meta-catalogue MCXC by 
Piffaretti et al. 2011). By contrast, although proposed about fif- 
teen years ago (e.g., Barbosa et al. 1996; Aghanim et al. 1997), 
it is only very recently that SZ- selected samples have reached 
a significant size, with publication of samples containing sev- 
eral hundred clusters from the Early SZ (ESZ) catalogue from 
the Planck Satellite (Planck Collaboration VIII 201 1), the South 
Pole Telescope (SPT, Reichardt et al. 2013) and the Atacama 
Cosmology Telescope (ACT, Hasselfield et al. 2013). 



^ Planck (http : //www . esa . int /Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



- Shopley Supercluster, (RA,DEC) = (202.6°,-31 .5°) 




Fig. 1 : The Shapley super-cluster as seen in the Planck survey. 
Upper panel: reconstructed thermal SZ map 3.2 x 1.8 square de- 
grees from MILCA (Hurler et al. 2010). The dotted circles rep- 
resent apertures of ^500 from the MCXC meta-catalogue around 
the resolved clusters. Lower panel: composite view of the opti- 
cal from DSS images. X-rays from ROSAT (light pink) survey 
and from the thermal SZ effect as seen in Planck (cyan). 



The usefulness of the different selection methods, particu- 
larly for cosmology, depends not just on the total number of clus- 
ters identified but also on how readily the selection function of 
the survey can be modelled, and on how well the observed clus- 
ter properties can be related to quantities such as the total cluster 
mass that are most readily predicted from theory (e.g., see Voit 
2005). It has proven difficult to capitalize on the large size of 
optical/IR cluster samples because the observable, the number 
of galaxies in each cluster, exhibits large scatter with respect to 
the total cluster mass (e.g., Johnston et al. 2007). In this regard 
the X-ray selected samples are considerably more powerful, due 
to the tighter correlations of X-ray properties with mass ( Amaud 
et al. 2005; VikhHnin et al. 2006; Pratt et al. 2009; Reichert 
et al. 2011; Maughan et al. 2012). Simulations predict that SZ- 
selected surveys may do even better, with a very tight relation 
between SZ signal and mass (e.g., da Silva et al. 2004; Motl 
et al. 2005; Nagai 2006; Wik et al. 2008; Aghanim et al. 2009; 
Angulo et al. 2012). Moreover, this relation, except at low red- 
shifts, corresponds to a nearly redshift-independent mass limit, 
thus allowing such surveys to reach to high redshift and provide 
a strong lever arm on growth of structure. 

We report on the construction and properties of the PSZ cat- 
alogue, which is to date the largest SZ-selected cluster catalogue 
and has value added through compilation of ancillary informa- 
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tion. It contains 1227 entries including many multiple systems, 
e.g., the Shapley super-cluster displayed in Fig. 1 together with 
a composite image. Of these 861 are confirmed, amongst which 
178 are new discoveries, whilst amongst the 366 candidate clus- 
ters 54 are of high reliability (CLASS 1 in our terminology), 170 
are reliable, and the remaining 142 are in the lowest reliability 
class. In Sect. 2 we start with a description of the Planck data 
used to provide cluster candidates, and the two different method- 
ologies (one of which has two independent implementations) 
used to carry out the extraction of the SZ sources. In Sect. 3 
we provide a characterization of the PSZ catalogue in terms of 
completeness, statistical reliability, and accuracy of cluster pa- 
rameters including size and photometry. Section 4 extensively 
describes validation of cluster candidates through pre-existing 
surveys and cluster catalogues in many wavebands, while Sect. 5 
describes the follow-up campaigns conducted by the Planck col- 
laboration to confirm new cluster discoveries. This leads to a 
description of the catalogue properties in Sect. 6. The physical 
properties of the clusters are exploited in Sect. 7. These include 
an update of the SZ-X-ray scaling relations from the Planck 
data, the measure of the X-ray flux for all SZ detections, and 
the production of homogenized SZ-mass estimates for 813 clus- 
ters with measured redshifts that are provided to the community 
as a value-added element to the Planck SZ catalogue. 

Throughout the article, the quantities M500 and R^m stand for 
the total mass and radius where the mean enclosed density is 500 
times the critical density at the cluster redshift. The SZ flux is de- 
noted 7500, where Y^m D\ is the spherically-integrated Compton 
parameter within 7?5oo, and Dp, is the angular-diameter distance 
to the cluster. The physical cluster quantities are computed with 
a fiducial ACDM cosmology with = 70 km s"^ Mpc"\ Dm = 
0.3 and fl^ = 0.7. Furthermore, all the fits are undertaken in the 
log-log plane using the BCES orthogonal regression method of 
Akritas & Bershady (1996), with bootstrap resampling, which 
allows for intrinsic scatter as well as uncertainties in both vari- 
ables. All uncertainties are given at 68 per cent confidence level 
and all dispersions are given in log 10. 

2. Construction of the Planck SZ Catalogue 

2. 1. Input Planck data 

The Planck catalogue of SZ sources is constructed from the total 
intensity data taken during the first 15.5 months of Planck survey 
observations. Raw data were first processed to produce cleaned 
time-lines (time-ordered information) and associated flags cor- 
recting for different systematic effects; channel maps were then 
produced for all the observing frequencies (see details in Planck 
Collaboration VI 2013; Planck CoUaboration II 2013). These 
maps, together with the associated beam characteristics, are the 
main inputs for the SZ-finder algorithms presented in Sect. 2.2. 
Following Planck Collaboration VIII (2011), we used the six 
highest- frequency Planck channel maps, from 100 to 857 GHz, 
to produce the catalogue of SZ detections. This optimizes the 
signal-to-noise (S/N) of the extracted SZ detections and the us- 
able sky fraction; see Appendix A for the choice of channel 
maps. 

In order to optimize the SZ detection, together with avoid- 
ing contamination of the PSZ catalogue by bright point sources 
(PS), the latter are masked from the channel maps prior to 
the SZ detection as detailed in the following. To construct the 
PS mask, we use the Planck Catalogue of Compact Sources 
(PCCS). The PCCS (Planck CoUaboration XXVIII 2013) is a 
collection of single-frequency source catalogues, one for each of 



the nine Planck frequency channels. The six single Planck-HFl 
frequency PS catalogues are used to first produce individual - 
frequency masks constructed by masking a radius equivalent to 
1.28 FWHM (3(Tbeam) around every point source detected with 
(S/N)ps > 10. Then a single common PS mask (see Fig. 2), 
which is the union of the six individual HFI-frequency channel 
masks, is constructed. It is applied to all six highest-frequency 
Planck channel-maps to mask the point sources prior to running 
the algorithms to detect SZ signal. The masked regions are filled 
using a harmonic in-painting method based on that of Bajkova 
(2005), which has the advantage of eliminating the discontinu- 
ities caused by the masking. In order to avoid any possible ar- 
tificial spurious detections at the edges of the in-painted area, 
we further reject detections within an expanded common mask, 
constructed using the same procedure as described above, but 
using a masking radius equivalent to 2.13 FWHM (5 (Xbeam) and 
covering less than 2.9% of the sky. 

Bright radio sources are known to exist at the centre of 
galaxy clusters, but they generally have steep spectra and hence 
their flux is significantly reduced at the six highest Planck fre- 
quencies where the PS mask is constructed and where the clus- 
ters are detected. The Perseus cluster (see Fig. 19 later and the 
associated discussion) is one exception, with a point source that 
is so bright that the cluster is masked and thus not included in 
the Planck SZ catalogue. 

2.2. Detection Methods 

The catalogue of SZ sources is the result of a blind multi- 
frequency search, i.e., no prior positional information on known 
clusters is used as input to the detection, by three detection algo- 
rithms briefly described below. These algorithms were described 
and tested using simulations (Melin et al. 2012). They were 
used to construct the Early SZ (ESZ) Planck sample by Planck 
Collaboration VIII (201 1). All three assume priors on the cluster 
spectral and spatial characteristics, which optimize the SZ de- 
tection by enhancing the SZ contrast over a set of observations 
containing contaminating signals. In the following we present 
the cluster model used as a template by the SZ-finder algorithms 
and we briefly describe the three detection methods (for details 
we refer the reader to Herranz et al. 2002; Melin et al. 2006; 
Carvalho et al. 2009, 201 1 ; Melin et al. 2012). 

2.2.1. Cluster model 

The baseline pressure profile model used in the detection meth- 
ods is the generalized NFW (Navarro et al. 1997) profile of 
Amaud et al. (2010). This profile model was constructed by 
combining the observed, scaled. X-ray pressure profile of 31 
clusters from the REXCESS sample (Bohringer et al. 2007) for 
R < Rsoo^, with the mean pressure profile from three sets of 
numerical simulations (Borgani et al. 2004; Nagai et al. 2007; 
Piffaretti & Valdamini 2008) for R500 < R < 57?5oo. New ob- 
servational constraints on the pressure distribution at R > R500 
have become available. Planck Collaboration Int. V (2013) con- 
strained the detection of the thermal pressure distribution out to 
about 3 7?5oo through stacking of the observed SZ profiles of 62 
nearby massive clusters detected with high significance in the 
Planck ESZ sample. The resulting profile is in agreement with 
that derived for the Coma cluster (Planck Collaboration Int. X 

^ Rsoo relates to the characteristic cluster scale R^ through the NFW 
concentration parameter C500 = 1.177 for the baseline profile (R^ = 

^500/C500)- 
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Fig. 2: The distribution, shown in Mollweide projection with the Galactic plane horizontal and the Milky Way centre in the middle, of 
the 1227 Planck clusters and candidates across the sky (red thick dots). The masked point-sources (black thin dots), the Magellanic 
clouds (large black areas) and the Galactic mask, covering a total of 16.3% of the sky and used by the SZ-finder algorithms to detect 
SZ sources, are also shown. 



2013). Both show a slightly flatter distribution in the outer parts 
(i.e., beyond 7?5oo) with respect to the predictions from the nu- 
merical simulations. These results are further confirmed by inde- 
pendent measurements from Bolocam in a smaller radial range 
(r < 27?5oo, Sayers et al. 2012b). Using the profile of Planck 
Collaboration Int. V (2013) does not affect the detection yield 
(see Sect. 3) and only slightly modifies the measure of the SZ 
flux density (see Sect. 7.5) as compared to the generalized NFW 
(GNFW) profile adopted in the three cluster. The fiducial model 
parameters for the GNFW profile are given by the parameteriza- 
tion of the pressure profile in Eq. 12 of Amaud et al. (2010). It 
states 



{csmx)y [1 + {csmxY]^-^^i- ' 

with the parameters 

[P^.csm.y.c^.p] = [8.40/z^q/^ 1.18,0.308, 1.05,5.49]. (2) 

The (weak) mass dependence of the profiles is neglected. Within 
the SZ-finder algorithms, the size and amplitude of the profile 
are allowed to vary but all other parameters are fixed. The cluster 
model is thus equivalent to a shape function characterized by 
two free parameters, its amplitude and a characteristic scale = 

^500 /^^500- 

2.2.2. Matched Multi-filter (MMF) 

Two different implementations of the matched multi-frequency 
filter algorithm (MMF 1 and MMF 3) are used to detect SZ clusters. 
Both are extensions, over the whole sky, of the MMF algorithm 



(Herranz et al. 2002; Melin et al. 2006). The matched filter op- 
timizes the cluster detection using a linear combination of maps 
(which requires an estimate of the statistics of the contamina- 
tion) and uses spatial filtering to suppress both foregrounds and 
noise (making use of the prior knowledge of the cluster pressure 
profile and thermal SZ spectrum). 

The MMF 1 algorithm divides the full- sky Planck frequency 
maps into 640 patches, each 14.66 x 14.66 square degrees, cov- 
ering 3.33 times the sky. The MMF 3 algorithm divides the maps 
into a smaller set of 504 overlapping square patches of area 
10 X 10 square degrees with the sky covered 1.22 times. The 
smaller redundancy of MMF 3 with respect to MMFl implies a 
potentially lower reliability of the SZ detections. In order to in- 
crease the reliability of the detections, the MMF 3 algorithm is 
thus run in two iterations. After a first detection of the SZ candi- 
dates, a subsequent run centred on the positions of the candidates 
refines the estimated S/N and candidate properties. If the S/N of 
a detection falls below the threshold at the second iteration, it 
is removed from the catalogue. For both implementations, the 
matched multi-frequency filter optimally combines the six fre- 
quencies of each patch. Auto- and cross-power spectra are di- 
rectly estimated from the data and are thus adapted to the local 
instrumental noise and astrophysical contamination, which con- 
stitutes the dominant noise contribution. Figure 3 illustrates, for 
a six arcmin filter size, the ensemble noise maps as measured 
by MMF 3 in each of the patches. For both MMFl and MMF 3, the 
detection of the SZ-candidates is performed on all the patches, 
and the resultant sub-catalogues are merged together to produce 
a single SZ-candidate catalogue per method. 
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Fig. 3 : The noise maps per detection patch of MMF 3 method mea- 
sured for a six arcmin filter. The noise ranges from 0.5 to 2 
times the average noise of the map, which is cry = 2.4 x 10""^ 
arcmin^. The Ecliptic polar regions, delimited by green con- 
tours, with increased redundancy in the observations define a 
deep survey zone covering in total 2.7% of the sky. It is less noisy 
than the areas near the Galactic plane, where the dust emission 
is higher. Two other zones are defined: a medium-deep survey 
zone of 41.3% coverage delimited by the red contours and with 
higher noise level; and a shallow- survey zone covering 56% of 
the sky and with the highest noise levels including regions near 
the Galactic plane. 



The candidate size in both algorithms is estimated by filter- 
ing the patches over the range of potential scales, and finding the 
scale that maximizes the S/N of the detected candidate. When 
merging the sub-catalogues produced from the analysis of indi- 
vidual patches, it is also the S/N of the detection (the refined S/N 
estimate for MMF 3) which is used when deciding which detec- 
tion of the candidate is kept. Furthermore, both MMF 1 and MMF 3 
can also be run with fixed cluster size and position to estimate 
the SZ signal. This version of the algorithms is used to assess 
the reliability of the association with known clusters and/or to 
refine the measurement of the integrated Compton parameters of 
known X-ray clusters, as presented in Sect. 7.2.1. 

2.2.3. PowellSnakes 

PowellSnakes (PwS) is different from the MMF methods. It is 
a fast Bayesian multi-frequency detection algorithm designed 
to identify and characterize compact objects buried in a dif- 
fuse background. The detection process is grounded in a statis- 
tical model comparison test. The statistical foundations of PwS 
are described in Carvalho et al. (2009), and more recently in 
Carvalho et al. (2011) with a greater focus on the Planck setup. 
PwS may be run either based on a Generalized Likelihood Ratio 
Test or in full Bayesian mode. This duality allows PwS mea- 
sured quantities to be consistently compared with those of the 
MMF algorithms. 

PwS also operates in square flat patches of 14.66 x 14.66 
square degrees. The total number of patches employed, of or- 
der 2800, varies with sky area but always guarantees a very 
large overlap; on average each cluster is detected about 4.7 
times. PwS detects putative clusters and at the same time it com- 
putes the evidence ratio and samples from the posterior distribu- 
tions of the cluster parameters. Then, it merges all intermediate 
sub-catalogues and applies the criterion of acceptance/rejection 



(Carvalho et al. 2011). PwS computes the cross-channel covari- 
ance matrix directly from the pixel data. To reduce the contam- 
ination of the background by the SZ signal itself, the estimation 
of the covariance matrix is performed iteratively. After an ini- 
tial estimate, all detections in the patch with S/N higher than 
the current target detection are subtracted from the data using 
their best-fit values and the cross-channel covariance matrix is 
re-estimated. This is PwS 'native' mode of background esti- 
mation that produces, on average, an S/N estimate about 20% 
higher than MMF. However, in order to produce a homogeneous 
Planck SZ catalogue from the three algorithms, it is possible to 
run PwS in 'compatibility' mode, skipping the re-estimation step 
to mimic more closely the evaluation of the background noise 
cross-power spectrum of the MMF algorithms and thus their eval- 
uation of the S/N. In this mode, PwS is a maximum likelihood 
estimator like the MMF. 

In the following, unless stated otherwise, all quoted or plot- 
ted S/N values from PwS are obtained in 'compatibility' mode 
in order to ensure homogeneity across the catalogue entries and 
in order to ease the comparison with the MMF outputs. 

2.3. Outputs of the detection methods 

Each of the three detection algorithms outputs a catalogue of SZ 
detections above S/N =4.5 outside the highest-emitting Galactic 
regions (this corresponds to a mask of about 15% of the sky, 
see masked area in Fig. 2) and the Small and Large Magellanic 
Clouds and outside the PS mask described in Sect. 2.1. The 
union PS-Galactic mask covers 16.3% of the sky. The survey 
area used for the SZ detections in Planck is thus 83.7% of 
the sky coverage. The three individual lists of SZ candidates 
are cleaned by removal of obvious false detections. These are 
spurious sources that pass the MMF and PwS filters despite the 
pre-processing step applied to the Planck channel maps, see 
Sect. 2.1. In order to identify them, we cross-match the SZ de- 
tections with an intermediate, low signal-to-noise cut of 4, cat- 
alogue of point sources detected at the highest frequencies of 
Planck. Galactic sources in dense and cold regions at high lat- 
itudes also contaminate the SZ detections outside the Galactic 
mask. These cold Galactic sources (CGS hereafter, see Planck 
Collaboration XXIII 201 1 ; Planck Collaboration XXII 201 1) are 
detected in the Planck channel maps following an optimized 
method proposed by Montier et al. (2010). The SZ detections 
matching with PS at both 545 and 857 GHz, or with CGS 
sources, all show a rising spectrum at high frequencies, indicat- 
ing that they are false detections. The SZ detections correspond- 
ing to such PCCS or CG sources are removed from the individual 
lists and from the published Planck catalogue of SZ sources. 

The three detection algorithms used in the present study de- 
ploy the GNFW cluster profile to detect SZ signal with the two 
parameters of the shape function, the central value and the char- 
acteristic scale let free, with = ^50o/<^50o- Each of the three 
algorithms therefore assigns, to each detected SZ candidate, a 
position with estimated uncertainty, a signal-to-noise value, and 
an estimated size, or equivalently ^500, with its uncertainty. 
The detection likelihood or the posterior probability of the inte- 
grated Compton parameter within 5^500, denoted ¥5^^^^^, exhibits 
a large correlation with the size. Figure 4 illustrates the like- 
lihood plots for two cases: a spatially-resolved cluster detected 
with a high signal-to-noise, Abell 2163; and a non-resolved clus- 
ter at high redshift (z - 1), PSZl G266.6-27.3 (also known as 
PLCK G266.6-27.3 in Planck CoUaboration XXVI 2011). We 
also show in Fig. 5 the distribution of maximum likelihood SZ 
fluxes (^sRsoo) and sizes (^500) for the MMF 3 detections. 
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Fig. 4: Illustration of the SZ size-flux degeneracy for two clusters detected by Planck. Right: Abell 2163 (S/N = 27) and left: PSZl 
G266.6-27.3 (S/N = 6 at z ^ 1). The contours show the 68, 95, and 99 percent confidence levels. 
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Fig. 5: Distribution of the maximum likelihood SZ flux Y^k^^^ 
and size ^500 for Planck SZ detections in the union catalogue 
down to S/N = 4.5. Detections associated with known or new 
confirmed clusters are shown as open black circles. SZ cluster 
candidates are shown as filled-red circles. 



This "degeneracy" between cluster size and SZ flux propa- 
gates the size uncertainty to the SZ flux estimate, increasing and 
biasing its value dramatically. This effect being so detrimental, 
both the SZ blind flux and size best-fit estimates, and respective 
error bars, are not quoted in the catalogue outputs to avoid their 
misuse. Only the full joint F5R5oo-^s, or equivalently ^srsoo-^soo, 



posterior probability contours provide a complete description of 
the information output by each detection method. They are thus 
provided for each detection. In order to use the flux measure, one 
ought to break the size-flux degeneracy. This can be achieved by 
a joint analysis with a high-resolution observation of the same 
objects, or by assuming a prior on, or fixing, the cluster size e.g., 
to the X-ray size. The SZ signal can then be re-extracted with 
an uncertainty much smaller than the variation of the joint Y-6 
probability distribution. 

We now perform a systematic comparison of the outputs 
of the three algorithms. We compare the S/N ratio. In addition 
and for purposes of illustration, we compare the best-fit blind 
Y value from maximum-likelihood or posterior probability out- 
puts, namely I^srsoo^- We show the comparison in Fig. 6, con- 
sidering detections down to S/N = 4.5. We quantify the differ- 
ence between a given quantity estimated by two different al- 
gorithms, Q2 and 2i, by fitting a power law to the data in the 
form 22/ep = 10^ {QilQ^Y with a pivot Gp = 6 for S/N and 
2p = 4 X 10"^ arcmin^ for Y^r^^q. The results are given in Table 1, 
including the scatter estimates. The raw scatter was estimated 
using the error- weighted vertical distances to the regression line. 
The intrinsic scatter on Fsoo was computed from the quadratic 
difference between the raw scatter and that expected from the 
statistical uncertainties. Table 1 also lists the mean difference in 
logarithm, A(log 2)=log(22/2i), computed taking into account 
both statistical errors and intrinsic scatter, estimated iteratively. 

2.3.1. Signal-to-noise 

A crucial ingredient of the SZ detection algorithms, either the 
MMFs or PwS, is the background cross-power spectrum used to 
estimate the noise level. It is evaluated from the data locally on a 
per-patch basis (see Fig. 3 for an example of the noise per patch 
across the sky). The algorithms, and implementations, slightly 

^ ^SRsoo rescaled to Ysqq for the fiducial GNFW model as 

75R500 = 1-79 X 7500 (Amaud et al. 2010). 
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Fig. 6: Comparison of S/N (top panels) and maximum likelihood Compton-parameter values (bottom panels) from the three detec- 
tion algorithms, MMFs and PwS, down to S/N = 4.5 after removing obvious false detections (see Sect. 2.3). In each panel, the red 
line denotes the equality line. The black line is the best fit to the data, and the dashed lines correspond to the ilcr dispersion about 
the fit relation. For clarity, error bars are omitted on Y^r^^q values in the plot, but are taken into account in the fit. The green line of 
slope fixed to unity corresponds to the mean offset between the two quantities. 



Table 1 : Quantitative comparison between S/N and maximum likelihood Compton-parameter Ysk^qq values from the three individual 
algorithms. Column(l): considered pair of algorithms ; Column(2-3): slope and normalization of the best-fit relation between the 
quantities estimated by the two algorithms, QilQ^ = (Qi/QpT, using BCES orthogonal regression, with the pivot being Qp = 6 
for S/N and 2p = 4 x 10"^ arcmin^ for ¥5^1^^^', Column (4-5) intrinsic and raw scatter around the best-fit relation; Column (6): mean 
difference in logarithm, A(log(2))=log(22/Gi); Column (7-8): corresponding intrinsic scatter and raw scatter. 
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MMF3-MMF1 
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+0.010 + 0.005 
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PwS-MMFl 


+0.041 + 0.004 


1.02 + 0.01 


0.04 + 0.01 


0.088 + 0.005 


+0.038 + 0.004 


0.040 + 0.007 


0.079 



differ with respect to the stabilization assumptions (e.g., smooth- 
ing) of the background noise cross-power spectrum and to the 
treatment of the background SZ signal, now acting as a contam- 
inant. These differences translate into variations in the S/N val- 
ues per method. In particular, when operated in ''compatibility'' 
mode (without background cluster subtraction), PwS estimation 
of the background cross-power spectrum is more affected than 
the MMF by SZ signal contamination. The SZ signal adds an 
extra component to the background noise producing lower S/N 
estimates. This is particularly noticeable when the SZ signal is 
very strong compared with background (typically S/N > 15). 



Despite the differences in background estimates, the yields 
from the three algorithms agree. In the left panel of Fig. 7, we 
show that the detection counts as a function of S/N for each de- 
tection method are in good overall agreement. The right panel 
of Fig. 7 shows the fraction of common detections over the 
union of detections from all three algorithms as a function of 
S/N. Sources with S/N > 8.5 are detected by all three meth- 
ods. However, we note that PwS number counts decrease more 
rapidly than MMF counts above S/N = 15. This reflects the be- 
haviour of PwS in "compatibility'' mode described above, which 
estimates a higher background than the MMF methods at high 
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Fig. 7: Left: detection number counts as a function of S/N of the individual algorithms. The S/N value in the union catalogue is that 
of the MMF3 detections when available, followed by that of PwS followed by MMFl (see Sect. 2.4). See text for discussion on the 
lower S/N values of PwS compared to the MMF -based algorithms. Right: fraction of common detection over counts from the union 
catalogue. Sources with S/N > 8.5 are detected by all methods. 



S/N. Figure 6 shows the comparison of the S/N estimates from 
all three methods. The agreement is good on average. The mean 
ratio (or the normalization at the pivot of the power-law rela- 
tion) deviates from unity by less than 2% and at less than 3 cr 
significance. Here again at high S/N values, we note the ten- 
dency for lower S/N in PwS as compared to MMF (Fig. 6), and 
indeed the slope of the power-law relation is smaller than unity 
{a = 0.94 + 0.01 for MMF 3). 

2.3.2. Photometry 

We now compare the best-fit Y values (from maximum like- 
lihood and posterior probability) for the three detection algo- 
rithms. The comparison (Fig. 6, lower panels) shows a system- 
atic bias with PwS, yielding slightly smaller values than MMF, 
typically by 10%. However, the slope is consistent with unity, 
showing that this bias is not flux dependent. The MMF values dif- 
fer from each other by less than 3% on average. The scatter be- 
tween Y estimates is dominated by the intrinsic scatter (Table 1). 
It is clearly related to the size-flux degeneracy, the ratio between 
Y estimates for a given candidate being correlated with the size 
estimate ratio, as illustrated by Fig. 8. The scatter becomes com- 
patible with the statistical scatter when a prior on the size is used, 
e.g., size fixed to the X-ray size. 

2.4. Definition of the Planck SZ catalogue 

As discussed above, the processing details of each algo- 
rithm/implementation differ in the computation of the back- 
ground noise. The significance of the detections in terms of S/N, 
although in overall agreement, differs from one algorithm to the 
other and translates into different yields for the candidate lists 
from the three algorithms. We choose to construct a catalogue 
of SZ candidates that ensures, through redundant detections, an 
increased reliability of the low S/N sources, when they are de- 
tected by two methods at least, together with maximizing the 
yield of the catalogue. 



The Planck SZ cluster catalogue described in the following 
is thus constructed from the union of the cleaned SZ-candidate 
lists produced at S/N > 4.5 by all three algorithms. It contains in 
total 1227 SZ detections above S/N = 4.5. Note that in order to 
ensure homogeneity, in terms of detection significance, the S/N 
values of PwS quoted in the the union catalogue are obtained in 
compatibility mode, whereas the S/N obtained from PwS native 
mode are quoted in the PwS individual list. The union catalogue 
is constructed by merging detections from the three methods 
within an angular separation of at most five arcmin, in agree- 
ment with Planck position accuracy shown later in Fig. 12. As 
mentioned, no reference photometry is provided. However a ref- 
erence position for the SZ detection is needed. For compatibil- 
ity with the ESZ Planck sample, in the case of matching de- 
tection between methods we arbitrarily choose to take the coor- 
dinates from the MMF 3 detection as the fiducial position (MMF 3 
was the reference method used to construct the ESZ Planck sam- 
ple). When no detection by MMF 3 above S/N = 4.5 is reported, 
we took the PwS coordinates as fiducial, and the MMFl coor- 
dinates elsewhere. The S/N values in the union catalogue are 
taken following the same order, which explains why the MMF 3 
curve in Fig. 6 coincides with the union curve. The cluster can- 
didates in the union catalogue are cross-referenced with the de- 
tections in the individual lists. The reference positions and the 
S/N values are reported in the union catalogue. Given the size- 
flux degeneracy, the full information on the degeneracy between 
size and flux is provided with each individual list in the form of 
the two-dimensional marginal probability distribution for each 
cluster candidate as discussed above. It is specified on a grid of 
256 X 256 values in 6s and Y^r^^q centred at the best-fit values 
found by each algorithm for each SZ detection. 



An extract of the Planck SZ catalogue is given in 
Appendix B. The full online table for union Planck catalogue, 
the individual lists of SZ detections, and the union mask used 
by the SZ-finder algorithms together with comments assembled 



8 



Planck Collaboration: Planck catalogue of Sunyaev-Zeldovich sources 

Table 2: Statistical characteristics of the Planck SZ catalogues. The intersection is defined as the set of SZ detections common 
between to three extraction algorithms. Completeness is a function of both 1^500 and ^500: the Ysoo at a given completeness is estimated 
by marginalizing over ^500, weighting each (^500, ^500) bin by the theoretically-expected cluster counts. Positional accuracy is the 
median angular separation between real and estimated positions. 



Stat. Reliability (%) ^500 10 ^ arcmin^ at C% completeness Positional accuracy (arcmin) 
C=50 C=80 C=95 



Union 84 0.61 1.2 3.2 1.2 

Intersection 98 0.85 1.8 6.6 1.1 

MMFl 87 0.75 1.6 4.7 1.2 

MMF3 91 0.71 1.5 3.8 1.2 

PwS 92 0.65 1.4 3.2 0.9 
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Fig. 8: Correlation between the ratio of i^sRgoo estimates with PwS 
and MMF3 and the ratio of size estimates, shown on a grid of 
sizes. 



in an external file are available at ESA's Planck Legacy Archive 
(PLA)4. 

3. Statistical Characterization 

The statistical characterization of the PSZ catalogue is achieved 
through a process of MonteCarlo quality assessment (MCQA) 
that can be applied to each individual catalogue and to the 
merged union catalogues. The statistical quantities produced in- 
clude completeness, fraction of detections associated with true 
clusters called, statistical reliability or purity, positional ac- 
curacy, and accuracy of parameter estimation. Together, these 
statistics describe the quality of detections in the catalogue. The 
quality of the parameter estimation, including astrometry (clus- 
ter position and extent), is determined through comparison with 

^ http : / /www .sciops.esa. int /index . php?page= 
Planck_Legacy_Archive&pro ject=planck. 



the parameters of the input clusters. The statistical characteris- 
tics of the different lists are summarized in Table 2. 

3. 1. MCQA Pipeline and simulations 

The MCQA pipeline contains a common segment producing 
simulated input catalogues and processed, source-injected maps, 
which are then fed into the detection pipeline. In summary, the 
pipeline steps per MonteCarlo loop are: 

1 . creation of an input cluster catalogue; 

2. injection of clusters into common simulated diffuse fre- 
quency maps, including beam convolution; 

3. injection of multi- frequency point sources; 

4. pre-processing of maps, including masking and filling point 
sources; 

5. detection and construction of individual cluster-candidate 
catalogues; 

6. construction of a union catalogue given merging criteria; 

7. collation of input and output catalogues, producing detec- 
tion truth-tables and catalogues of unmatched spurious de- 
tections^. 

To estimate the completeness, clusters are injected into the 
real data. In this case, steps 3 and 4 are skipped and each detec- 
tion algorithm estimates noise statistics on the real data prior to 
injection in order to avoid artificially raising the S/N and biasing 
the completeness estimates. The pressure profiles of the injected 
clusters follow that described in Sect. 2.2.1. To account for the 
profile variation across the cluster population, the profile param- 
eters are drawn from the covariance matrix of the 62 measured 
pressure profiles from Planck Collaboration Int. V (2013), ensur- 
ing that the injected profiles are consistent with measured disper- 
sion and consistent, on average, with the extraction filter. The in- 
jected clusters are convolved with effective beams in each pixel 
including asymmetry computed following Mitra et al. (201 1). 

The simulated input cluster catalogues differ for statistical 
reliability and completeness determination. For completeness, 
clusters injected in real data are drawn from a uniform distribu- 
tion in (Ysoo, 6500) so as to provide equal statistics in each com- 
pleteness bin. To avoid an over-contamination of the signal, in- 
jected clusters are constrained to lie outside an exclusion radius 
of 57?5oo around a cluster, either detected in the data or injected. 

For the statistical reliability estimation of the input cluster 
distribution injected in simulations is such that cluster masses 
and redshifts are drawn from a Tinker et al. (2008) mass func- 
tion and converted into the observable parameters (Isoc^soo) 

^ A cluster is considered to be matched if there is a detection within 
five arcmin of its position. 
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using the Planck ESZ Ysoo-Msoo scaling relation (Planck 
Collaboration X 201 1). The simulated maps consist of CMB re- 
alizations, diffuse Galactic components and instrumental noise 
realizations, including realistic power spectra and inter-detector 
correlations, from the FFP6 simulations (Planck Collaboration 
XII 2013; Planck Collaboration ES 2013). Residual extragalactic 
point sources are included by injecting, mock-detecting, mask- 
ing and filling realistic multi-frequency point sources using the 
same process as for the real data (see Sect. 2. 1). 

3.2. Completeness 

The completeness is the probability that a cluster with given 
intrinsic parameters (^500,^500) is detected given a selection 
threshold (here in S/N). 

If the Compton-F estimates are subject to Gaussian errors, 
the probability of detection per cluster follows the error func- 
tion and is parameterized by (Tn(^5oo), the standard deviation of 
pixels in the multi-frequency matched- filtered maps for a given 
patch / at the scale ^500, the intrinsic Compton and the de- 
tection threshold q: 



P(d\Y5oo,(TYi(e5oo),q)=^ 



1 -h erf 



^500 ~ Qcr YijOsoo) 



(3) 



where erf(x) = (2/;:) J^^ exp (-^^) dt and d is the Boolean detec- 
tion state. 

The completeness of the catalogue, thresholded at S/N q, is 
expected to follow the integrated per-patch error function com- 
pleteness 



C(y500, ^500) - ^ /sky,/^(<^l^^500, CT r/(^50o), q) , 



(4) 



where /sky,/ is the fraction of the unmasked sky in the patch /. The 
true completeness departs from this theoretical limit. This is due 
to the non-Gaussian nature of the noise dominated by the astro- 
physical, namely Galactic, contamination. This is also the case 
when the actual cluster pressure profile deviates from the GNFW 
used in the SZ-finder algorithms, or when the effective beams 
deviate from constant symmetric Gaussians, and also when the 
detection algorithm includes extra steps of rejection of spurious 
sources not formulated in Eq. 3. This is why an MCQA-based 
assessment of the completeness is essential to characterize the 
Planck detections. 

The MonteCarlo completeness of each of the individual lists 
and the union catalogue are shown in Fig. 9. The MMF lists are 
consistent with one another at ^500 > 4 arcmin, but MMF 3 is 
more complete at lower radii. This is due to an extra step imple- 
mented in MMF 1 that rejects as spurious the detections estimated 
to be point-like. The union improves upon the completeness of 
each of the individual catalogues, because it includes the faint 
real detections by one method alone. In contrast, the intersec- 
tion of the lists from the three algorithms, while more robust, 
is markedly less complete than the union and each of the in- 
dividual catalogues. The intersection and union catalogues rep- 
resent the extremes of the trade-off between statistical reliability 
and completeness. The quantities for each of the catalogues, plus 
the union and intersection, are summarized in Table 2. Figure 9 
shows four constant ^500 slices through the completeness con- 
tours for MMF 3 , comparing the MCQA-based completeness with 
the integrated error function completeness. At radii smaller than 
6 arcmin, the MCQA-based completeness is systematically less 
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Fig. 9: Top panel: differential completeness as a function of 
(5^500,^500) for each detection algorithm (MMFl in blue, MMF 3 
in red, and PwS in green) and for the union (shaded area) and 
intersection (black) catalogues. From bottom to top, the solid, 
dashed, and dotted lines show 15%, 50% and 85% complete- 
ness, respectively. Bottom panel: slices through the MCQA- 
based completeness function at various ^500 for MMF 3 compared 
to the error function approximation (solid curves). 



complete, and the drop-off of the completeness function shal- 
lower, than the theoretical expectation. This effect is a con- 
sequence of the variation of intrinsic cluster profiles from the 
GNFW profile assumed for extraction. 



10 



Planck Collaboration: Planck catalogue of Sunyaev-Zeldovich sources 



100 - 



- 95 



.CO 

DC 
> 



O 85 



90 - 




Union 
Intersection 
IVIIVIF1 
IVIIVIF3 
PwS 



80 I I I I I I I I I I I I I I I I I I I I ■_ 



100 - 



4 6 8 10 12 14 

S/N 

Fig. 10: Cumulative statistical reliability, defined as the fraction 
of sources above a given S/N associated with a "real" cluster 
from the simulated input catalogue. 



3.3. Statistical reliability 

The fraction of detections above a given S/N that are associ- 
ated with a real cluster is characterized by injecting clusters into 
high-fidelity simulations of the Planck channels. Unassociated 
detections from these simulations define the fraction of spuri- 
ous detections. We have verified that the simulations produced 
detection noise cty^q^ consistent with the real data and that the 
simulated detection counts match the real data. 

The cumulative fraction of true clusters, as characterized by 
the simulations, is shown for the output of each detection algo- 
rithm and for the union catalogue in Fig. 10. The union catalogue 
is less pure than any of the individual lists because it includes all 
the lower-reliability, individual-list detections, in addition to the 
more robust detections made by all three SZ-finder algorithms. 
The union catalogue constructed over 83.7% of the sky at S/N of 
4.5 is 84% pure. 

The fraction of false detections is dominated by systematic 
foreground signals, in particular Galactic dust emission. This is 
illustrated in Fig. 1 1 by the effect of dust contamination on the 
cumulative reliability. We define two sky regions by the level 
of dust contamination: "region 1" is the low dust-contamination 
region outside of the Planck Galactic dust, and PS, mask that ex- 
cludes 35% of the sky. This mask is used in Planck Collaboration 
XX (2013) for cosmological analysis of SZ counts. "Region 2" 
is the complementary region included by the smaller 15% dust 
mask but excluded by the 35% mask. When the larger Galactic 
dust mask is applied leaving 65% of the Planck sky survey in 
which to detect SZ signal, the statistical reliability increases 
from 84% in 83.7% of the sky to 88% in 65% of the sky. As 
seen in Fig. 1 1 upper panel, the reliability of the detections dete- 
riorates markedly in "region 2" relative to "region 1". The noisy 
behaviour of the curves in Fig. 1 1 upper panel is due to the re- 
duced size of sky area used in the analysis. 
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Fig. 1 1 : Top panel: cumulative reliability for the union and inter- 
section catalogues, as a function of dust contamination. Region 
1 is the low-dust contamination region, being the 65% of the sky 
outside the Galactic dust mask, and region 2 is the complemen- 
tary dustier region added to this when the smaller 15% dust mask 
is applied. The Gaussian noise limit is the expected reliability 
from purely Gaussian fluctuations. Bottom panel: histogram of 
the J- signal in a typical filtered patch from a null-test simula- 
tion, compared to the best-fit Gaussian (black dashed line). The 
distribution of j-noise is non-Gaussian. 



In both regions, the spurious count much higher than is pre- 
dicted by Gaussian fluctuations. This reflects the non-Gaussian 
nature of the filtered patches. The bottom panel of Fig. 1 1 il- 
lustrates this for a typical mid-latitude patch from a null-test 
simulation with no injected clusters. The patches are well ap- 
proximated as Gaussian at deviations smaller than 3 cr (consis- 
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Fig. 12: Distributions of positional error for each catalogue, nor- 
malized by the total number of detections in the catalogue. By 
construction, the positional error is defined to be less than five 
arcmin. 



tent with the assumptions of Eq. 3), but show enhanced numbers 
of high significance deviations, which can translate into spurious 
detections. 

3.4. Positional Accuracy 

Positional accuracy is characterized by the radial offset between 
estimated and injected positions. The distribution of position er- 
ror is shown in Fig. 12, for each individual list and the union cat- 
alogue. In contrast to the MMFs, which estimate the maximum- 
likelihood position, the PwS position estimator is the mean of 
the position posterior, which produces more accurate positional 
constraints. The union catalogue positions are taken from MMF3 
if available, followed by PwS and then MMFl. Its positional esti- 
mates are hence consistent with the MMFs. The mode of the union 
distribution is consistent with a characteristic position error scale 
of half an HFI map pixel (0.86 arcmin). 
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Fig. 13: Distributions of the ratio of detected over injected pa- 
rameters for F5R500 and ^500. 



3.5. Parameter Recovery 

The Compton Y^r^^q is characterized by comparing detected and 
input values for matched detections from the injection of clusters 
into the real data (see Fig. 13). The injection follows the scheme 
outlined above with one exception: input cluster parameters are 
drawn using the Tinker mass function and the scaling relations 
discussed above for reliability simulations. This ensures a real- 
istic distribution of parameters and S/N values. 

What we characterize is slightly different for each catalogue. 
For the MMFs, we characterize the maximum-likelihood point of 
the 2-D degeneracy contours provided in the individual lists. For 
PwS, we characterize the mean of the marginal distribution for 
each parameter. In each case, the 2D (I^srsoq, ^s) are marginalized 
over position. The contours are scaled for each cluster and are 
time consuming to compute, so we characterize the parameters 



from a lower-resolution grid that is better suited to Monte- Carlo 
analysis.^ 

The scatter between input and detected parameters is shown 
in Fig. 14 as an example for PwS. Biases are evident at both the 
low and high end for Y^r^^q . The low-flux bias is the Malmquist 
bias related to the S/N > 4.5 threshold. The high-flux bias is due 
to a hard prior on the upper limit for cluster radius. Figure 14 
also shows the distribution of the ratio of estimated over injected 
parameters. The median and median absolute deviation of these 
ratios are shown in Table 3. 



^ PwS does not resort to a low-resolution scale grid and always works 
at the full resolution. 



12 



Planck Collaboration: Planck catalogue of Sunyaev-Zeldovich sources 




lO"" 10'^ 10'" 10"' 10° 10° 10^ 10^ 

Y^pj Injected [arcmin^] e^^^ Injected [arcmin] 

Fig. 14: Injected versus detected values of Y^r^^^q (left panel) and ^500 (right-panel), illustrated for PwS. 



Table 3: Median and median absolute deviation (MAD) of the 
ratio of detected over injected parameters. 





^5R500 

median 


MAD 


e 

median 


MAD 


MMFl .... 


1.09 


0.39 


1.17 


0.70 


MMF3 


1.02 


0.34 


1.19 


0.69 


PwS 


0.99 


0.27 


1.21 


0.56 



The distributions for flux are positively skewed due to 
Malmquist bias. The median ratios of the flux recoveries are con- 
sistent with unity for MMF3 and PwS and are slightly higher for 
MMF 1 . The recovery of is biased high in the median by about 
20% for each of the codes. This bias is a consequence of the in- 
trinsic cluster profile variation and disappears when the injected 
profiles match the detection filter. The Ysr^^^ estimate by con- 
trast is relatively unaffected by profile variation. The parameter 
constraints from PwS are tighter than the MMFs due to the PwS 
priors and the definition of the estimator as the expected value 
of the parameters rather than the maximum likelihood. 

4. External Validation 

The cluster-candidate catalogue constructed from the union of 
all three SZ-finder algorithms undergoes a thorough validation 
process that permits us to identify previously-known clusters and 
to assess the reliability of the Planck SZ candidates not associ- 
ated with known clusters. In order to achieve this, we make use 
of the existing cluster catalogues and we also search in optical, 
IR, and X-ray surveys for counter-parts at the position of the 
Planck SZ sources. In Sect. 5, we present the follow-up pro- 
grammes that were undertaken by the Planck collaboration in 
order to confirm and measure the redshifts of the Planck candi- 
date new-clusters. 



The first step of the validation of the PSZ catalogue is to 
identify among the Planck SZ candidates those associated with 
known clusters. For this purpose, we use existing X-ray, opti- 
cal or SZ cluster catalogues. A positional matching is not suffi- 
cient to decide on the association of a Planck SZ source with a 
previously-known cluster, and a consolidation of the association 
is needed. For the X-ray associations, a mass proxy can be built 
and used to estimate the SZ flux, S/N, etc, that are compared with 
measured quantities for the Planck cluster candidates. In con- 
trast to the X-ray clusters, optical clusters either have no reliable 
mass estimates or suffer from large uncertainties in the mass- 
richness relations. In this case, the consolidation cannot be per- 
formed uniquely through the coherence of measured versus pre- 
dicted properties. It rather relies on extra information from sur- 
veys in the X-ray, optical, or IR at the Planck cluster-candidate 
positions. 

In the following, we detail the search for counter-parts in op- 
tical, IR, and X-ray surveys; list the cluster catalogues used for 
the identification; and finally present the identification procedure 
followed to associate Planck SZ detections with bona fide clus- 
ters. In this process, we define quality flags for the association of 
Planck SZ detections with external data. We set 2 = 1 for high- 
reliability associations, i.e., very clear cluster signatures, 2 = 2 
for reliable associations, and 2 = 3 for low-reliability associa- 
tions, i.e., unclear cluster signature. 

4. 1. Search for counter-parts of Planck detections in 
surveys 

We made use of the ROSAT All Sky Survey (RASS, Voges et al. 
1999), the all-sky survey with the Wide-field Infrared Survey 
Explorer (WISE, Wright et al. 2010), and the Sloan Digital Sky 
Survey (SDSS, York et al. 2000) to search for counter-parts of 
the Planck SZ detections. This information was used in two 
ways. When Planck detections were associated with known clus- 
ters from catalogues, in particular in the optical, the counter- 
parts in RASS, WISE, or SDSS helped in consolidating the asso- 
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elation, increasing the confidence in the identification of Planck 
candidates with known clusters. When no association between 
Planck detections and previously-known clusters was found, the 
information on the counter-parts, in the surveys, of Planck SZ 
detections was used to assess the reliability of the Planck cluster 
candidates, i.e., clear or unclear cluster signatures. 

4.1.1. Search in RASS data 

As detailed in Planck Collaboration Int. IV (2013), the vaUda- 
tion follow-up with XMM-Newton has shown the importance of 
the RASS data to assess the reliability of the Planck sources. 
In particular, Planck Collaboration Int. IV (2013) showed that a 
large fraction of Planck clusters are detectable in RASS maps, 
but this depends on the region of the sky and on the ratio Fsoo/^x 
which exhibits a large scatter (see later in Fig. 31 the case of the 
PSZ sources). We therefore exploit the RASS data to consoli- 
date the identification with clusters from optical catalogues (see 
below Sect. 4.3.2) and to assess the reliability of the Planck SZ 
candidates. 

We first perform a cross-match with the RASS bright source 
catalogue (BSC, Voges et al. 1999) and the faint source cata- 
logue (FSC, Voges et al. 2000) within a five-arcmin radius of 
the position of each of the Planck SZ detections. We then per- 
form a reanalysis of the RASS data following the methodology 
and prescriptions given by Bohringer et al. (2000, 2004) and 
Reiprich & Bohringer (2002). We compute count-rate growth 
curves in order to check for the extension of the signal. We es- 
timate the source flux from both the growth curve (when ade- 
quate) and from a fixed five-arcmin aperture radius with respect 
to the surrounding background (after PS subtraction). We then 
derive the associated signal-to-noise in RASS, (S/N)rass- For 
this, we make use of the RASS hard-band, [0.5-2] keV, data that 
maximize the S/N of the detections. We furthermore computed 
the source density map of the BSC and FSC catalogues and the 
associated probability that a Planck cluster candidate will be as- 
sociated with a B/FSC source within a radius of five arcmin. For 
the BSC, the probability of chance association is relatively low, 
with a median <1 %. As detailed in Planck Collaboration Int. 
IV (2013), the correspondence of a Planck SZ-candidate with 
a RASS-BSC source is a semi-certain association with a real 
cluster, whereas for the FSC catalogue the probability of chance 
association is larger, 5.2%. 

We define a quality flag, Grass, for the association of Planck 
candidates with RASS counter-parts using both the signal-to- 
noise in RASS and the association with B/FSC sources. This is 
of particular importance for the Planck candidate new clusters. 
Based on the results from Planck Collaboration Int. IV (2013), 
the quality of the association with RASS counter-parts is high. 
Grass = 1, for Planck cluster candidates matching a RASS- 
BSC source or with (S/N)rass > 2. We find a total of 887 out of 
1227 Planck SZ detections in this category, with mean and me- 
dian signal-to-noise of 7.4 and 5.8, respectively. The quality is 
poor. Grass = 3, for RASS counter-parts with (S/N)rass < 0-5 
in regions of reasonable depth (quantified by the probability of 
chance association with FSC sources being larger than 2.5% 
(Planck CoUaboration Int. IV 2013)). 

4.1.2. Search in SDSS data 

We performed a systematic search for counter-parts in the SDSS 
Data Release DR9 (SDSS-III Collaboration et al. 2012) at the 
position of all the Planck SZ detections. This was performed 



based on a cluster-finder algorithm developed by Fromenteau 
et al. (2013) to search for red galaxy over-densities in the SDSS 
galaxy catalogues. 

For each associated counter-part within a five arcmin cir- 
cle centred at the position of the Planck SZ detection, a quality 
criterion is defined on the basis of a fit to the luminosity func- 
tion and the associated mass limit, and on the number of galax- 
ies within five arcmin, A^gai, such that we have GsDSS,dat = 1^ 
i.e., high quality, for cases where A^gai > 40 and for masses 
M200 > 5.7 X 10^"^ Mq, GsDSS,dat = 2, i.e., good quality, for 
A^gai between 40 and 20 for masses between 1.5 x 10^"^ M© and 
5.7 X 10^"^ Mo, and GsDSS,dat = 3 otherwise. 

The cluster-finder algorithm outputs the position of the 
counter-part (Brightest Cluster Galaxy (BCG) and bary centre) 
and the estimated photometric redshift. When spectroscopic data 
are available for the brightest selected galaxy a spectroscopic 
redshift is also reported. The outputs of the cluster-finder algo- 
rithm are compared to those obtained by Li & White (2013) from 
different method based on the analysis of the full photometric- 
redshift probability distribution function (Cunha et al. 2009). In 
this approach, the position and redshift in the SDSS data that 
maximizes the S/N are considered as the best estimates for the 
counter-parts of the Planck SZ detections. 

4.1.3. Search in WISE data 

The Wide-field Infrared Survey Explorer (WISE, Wright et al. 
2010) provides an all-sky survey at 3.4, 4.6, 12, and 22 jjm (Wl, 
W2, W3, W4) with an angular resolution of 6.1 to 12.0 arcsec in 
the four bands. 

We search for counter-parts of the Planck SZ detections in 
the WISE source catalogue in two ways. On the one hand, we 
run an adaptive matched filter cluster finder developed by Aussel 
et al. (2013), similar to the one described by Kepner et al. (1999), 
using the cluster members' luminosity function of Lin et al. 
(2012). The background counts were determined from the neigh- 
bouring square degree in the vicinity of the Planck cluster can- 
didate, excluding regions of fifteen arcmin centred on candidate 
positions. On the other hand, we use a method developed by 
Aghanim & Fromenteau (2013) based on a search for overdensi- 
ties of bright (Wl < 17) and red (W1-W2 > 0) sources within a 
five-arcmin radius circle centred on the position of Planck detec- 
tions with respect to a background computed in a fifteen-arcmin 
radius area. 

Aghanim & Fromenteau (2013) find that a good-quality as- 
sociation between a Planck SZ-detection and a counter-part 
overdensity in WISE data is reached when there are at least ten 
galaxies above 2 cr in the five-arcmin search region, and when 
the corresponding fraction of galaxies is at least 30% of the 
total number of galaxies retained in the fifteen-arcmin circle. 
Performing the search for counter-parts of an ensemble of ran- 
dom positions on the sky, we compute the purity of the detec- 
tions, i.e., the probability of a Planck candidate having a real 
counter-part in the WISE data as opposed to a chance associ- 
ation. The quality criterion for the association between Planck 
detection and WISE overdensity is high, Gwise = 1 , for a pu- 
rity larger than 90%. When it Hes between 90% and 80% the 
association of Planck SZ-detections and WISE overdensities is 
assigned a lower quality criterion Gwise = 2. We set the quality 
of the association to Gwise = 3, bad, when the purity is be- 
low 80%. We find 856 Planck SZ detections with high or good 
quality counter-parts in WISE data, including 658 Gwise = 1 
detections. 
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4.1.4. DSS images 

Finally for each Planck cluster candidate, the second Digitized 
Sky Survey^ (DSS) database was queried for a field of 5 x 

5 arcmin^ centred at the position of the Planck SZ detections 
in the r and ir bands. The DSS images were used for visual in- 
spection.^ Clusters and rich groups out to z - 0.3 to 0.4 can 
easily be identified in these plates as an obvious concentration 
of galaxies. This qualitative information was thus used: (i) to 
consolidate some identifications of Planck SZ detections with 
previously-known clusters; (ii) to optimize our strategy for the 
follow-up observations of Planck candidates (see Sect. 5); and 
(iii) to qualitatively assess the reliability or significance of the 
Planck SZ detections. 

4.2. Cluster catalogues 

We now present the ensemble of catalogues that were used to 
identify the Planck SZ detections with previously-known clus- 
ters. In the case of the ROSAT- and SDSS-based catalogues, we 
have used homogenized quantities, see below, that allowed us to 
perform the identification with comparable association criteria, 
which ensures homogeneity in the output results. 

MCXC meta-catalogue - For the association of Planck SZ 
candidates with previously-known X-ray clusters, we use the 
Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC, 
Piffaretti et al. 201 1 , and reference therein) constructed from the 
publicly available ROSAT All Sky Survey-based and serendip- 
itous cluster catalogues, as well as the Einstein Medium 
Sensitivity Survey. For each cluster in the MCXC several prop- 
erties are available, including the X-ray coordinates, redshift, 
identifiers, and standardized luminosity, Lx,5oo, measured within 
RsoQ. The MCXC compilation includes only clusters with avail- 
able redshift information (thus X-ray luminosity) in the origi- 
nal catalogues. We updated the MCXC, considering the first re- 
lease of the REFLEX-II survey (Chon & Bohringer 2012), the 
third public release of clusters from the MACS sample (Mann 

6 Ebeling 2012), individual MACS cluster publications and a 
systematic search in NED and S 1MB AD for spectroscopic red- 
shift for clusters without this information in the ROSAT cata- 
logues. This yields an ensemble of 1789 clusters with z and 
^x,5oo values, adding 20 MACS clusters, 21 REFLEX-II clus- 
ters and 5 SGP clusters to the MCXC. For these clusters, the 
expected Compton-parameter, Y^^^, and size, 6^^^, are estimated 
combining the M5oo-^x,50o relation of Pratt et al. (2009) and the 
^500-^500 relation given by Arnaud et al. (2010). The expected 
S/N ratio, (S/N)^^, is computed taking into account the noise 
within ^^QQ at the cluster location. We furthermore supplement 
the updated MCXC with 74 clusters from ROSAT catalogues 
without redshift information and 43 unpublished MACS clusters 
observed by XMM-Newton or Chandra. For these 117 objects, 
only centroid positions are available. Finally, we considered the 
published catalogues from XMM-Newton serendipitous cluster 
surveys with available redshifts, the XCS catalogue (Mehrtens 
et al. 2012), the 2XMMi/SDSS catalogue (Takey et al. 201 1) and 
the XDCP catalogue (Fassbender et al. 2011). However, these 
catalogues mostly extend the MCXC to lower masses and only 
two Planck candidates were found to be associated with these 
new clusters. 

http://stdatu.stsci.edu/dss/. 
^ Images from the RASS, SDSS and WISE surveys at the position of 
the Planck SZ detections were also inspected. 




10 100 

Fig. 15: Mass to richness scaling relation, Msoq-^l*, for the 444 
MCXC clusters included in the WHL12 catalogue (Wen et al. 
2012). The best-fit relation, from BCES fit, is given by the solid 
blue line. We adopted 15% uncertainties on the MCXC masses 
as prescribed in Piffaretti et al. (201 1). As no uncertainty is pro- 
vided for the WHL12's richness, we arbitrarily assumed a 20% 
uncertainty for all richness values. The blue shaded area shows 
the associated errors on the best-fit, while the dashed line marks 
the intrinsic scatter. 



Optical-cluster catalogues - The identification of the Planck 
SZ candidates with clusters known in the optical is based on the 
Abell (Abell 1958) and the Zwicky (Zwicky et al. 1961) cluster 
catalogues. Furthermore, we have used four different catalogues 
of clusters based on the Sloan Digital Sky Survey (SDSS, York 
et al. 2000) data: (1) the MaxBCG catalogue (13,823 objects, 
Koester et al. 2007); (2) the GMBCG catalogue (55, 424 objects, 
Hao et al. 2010); (3) the AMF catalogue (69, 173 objects, Szabo 
et al. 2011); and (4) the WHL12 catalogue (132,684 objects. 
Wen et al. 2012). We refer the reader to Wen et al. (2012) for 
a comparison of the existing SDSS-based catalogues of clusters 
and groups. Each of the SDSS-based catalogues provides an es- 
timated richness; we first start by homogenizing the richness es- 
timates to that of WHL12. For each catalogue, we compute the 
median ratio of WHL12's richness to that of the considered cat- 
alogue over its intersection with WHL12's. We then renormalize 
the individual richness by the corresponding ratio. The correct- 
ing factors applied to the richness estimators^ are respectively 
1.52, 1.75, and 0.74 for MaxBCG, GMBCG, and AMF, obtained 
from 7627, 17245, and 1358 common clusters. The richness is 
then related to the halo mass, M500, by extending the Wen et al. 
(2012) richness-mass relation provided on about 40 clusters to 

9 Field NGALS_R200 for MaxBCG, GM_SCALED_NGALS for 
GMBCG and LAM200 for AMR 

We considered the associations of clusters with positions match- 
ing within 6 arcsec radius and with Az < 0.05 (typical uncertainty for 
photometric redshifts in SDSS). 

^ ^ Their M200 are taken from the literature either from weak lensing or 
X-ray measurements (Wen et al. 2010). 
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Fig. 16: Identification of the Planck cluster candidates with X-ray clusters from the MCXC catalogue. Black points are candidates 
firmly identified with MCXC clusters, while green points are candidates with no association. Left panel: distance of the Planck 
position to the position of the closest MCXC cluster as a function of the distance normalized to the cluster size ^^q. Middle panel: 
S/N normalized to the expected value as a function of normalized distance. Right panel: SZ flux, y5oo,psx, re-extracted fixing the 
position and size to the X-ray value, as a function of expected values. The red line is the equality line. In all panels, 7^^^, and 6^^^ 
are estimated from the cluster X-ray luminosity used as mass proxy (see text). 



444 MCXC clusters, with masses estimated from the X-ray lumi- 
nosities. The data points and the best-fit scaling relation are pre- 
sented in Fig. 15. The derived Msqq-Ri* and Lx,200-^l* relations 
are compatible with the findings of Wen et al. (2012). We find 
log(M5oo/lO^'^M0) = (-2.00±0.17)-h(1.37±0.10)xlog7?L*.The 
relation presents a large intrinsic log- scatter, (Tint = 0.27 + 0.02, 
hampering any accurate estimation of the cluster mass. This is 
further illustrated by the richest clusters with > HO hav- 
ing MCXC masses systematically below the best-fit M2oo-^l* 
relation (although within the Icr intrinsic scatter). 

SZ catalogues - At millimetre wavelengths, we cross-check 
the Planck SZ catalogue with the recent ACT and SPT samples 
(Menanteau et al. 2010; Vanderlinde et al. 2010; WilHamson 
et al. 2011), including the most recent data that increased the 
number of SZ detections and updated the redshift estimates for 
the clusters (Reichardt et al. 2013; Hasselfield et al. 2013). We 
have furthermore identified the Planck SZ detections associated 
with previous SZ observations of galaxy clusters from the lit- 
erature. We used a compilation of SZ observations conducted 
with the numerous experiments developed during the last 30 
years (Ryle, OVRO, BIMA, MITO, Nobeyama, SZA, APEX- 
SZ, AMI, Diabolo, Suzie, Ryle, AMIBA, ACBAR, etc.). 

4.3. Identification with previously-known clusters 

4.3.1. Identification with X-ray clusters 

The Planck SZ candidates are cross-checked against previously- 
known X-ray clusters from the updated version of the MCXC. 
For a given Planck candidate-cluster we identify the closest 
MCXC cluster. The reliability of the association is assessed 
based on distance, D, compared to the cluster size and on the 
measured Ysoo and S/N values compared with the expected val- 
ues (see Fig. 16). Two clouds of points stand out in the scat- 
ter plot of absolute versus relative distance, ^/^qo ^^^S- 16, left 
panel). They correspond to two clouds in the scatter plot of the 

The information of the second closest is also kept to identify poten- 
tial confusion or duplicate associations. 



measured over expected S/N ratio versus D/^^JFig. 16, middle 
panel). 

The association process follows three main steps. First, we 
provisionally assign an X-ray identification flag based on dis- 
tance: 

- Qx = 3ifD > le^l^ and Z) > 10 arcmin. Those are consid- 
ered as definitively not associated with an MCXC cluster in 
view of Planck positional accuracy and cluster extent. 

arcmin. Those are associ- 
ated with an MCXC cluster. 

- 2x = 2 otherwise, corresponding to uncertain associations. 

We then refine the classification. In the 2x = 1 category, we 
identify outliers in terms of the ratio of measured to expected 
S/N and Fsoo, taking into account the scatter and the size-flux 
degeneracy. Their flags are changed to 2x = 2. In some cases, 
two distinct 2x > 1 candidates are associated with the same 
MCXC cluster. The lowest S/N detection is flagged as 2x = 2. 

In the final step, we consolidate the status of 2x < 3 can- 
didates. We first re-extract the SZ signal at the X-ray position, 
both leaving the size free and fixing it at the X-ray value. The 
7500 obtained with the cluster and size fixed to the X-ray values 
are compared to the expected values, 7^^^, in the right panel of 
Fig. 16. For bona fide association, we expect no major change 
of Fsoo and S/N, with, on average, a better agreement with the 
expected Fsoo value and some decrease of S/N. 

- For 2x = 1 candidates, the re-extracted Fsoo and S/N values 
are compared to both blind and expected values (as a func- 
tion of distance, S/N, etc.) to identify potential problematic 
cases, e.g., important decrease of S/N or outliers in terms of 
measured-over-expected Y^m ratio. We found only one such 
case, whose flag is changed to 2x = 2. The identification of 
other candidates is considered as consolidated, with defini- 
tive flag 2x = 1. 

- We then examine the 2x = 2 candidates. We consider the 
re-extracted Y^m and S/N, but also perform a visual inspec- 
tion of the SZ maps and spectra and ancillary data, including 
RASS and DSS images. The 2x = 2 candidates were iden- 
tified as clearly identified as multiple detections of extended 
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clusters or duplicate detections of the same clusters by differ- 
ent methods that were not merged (the former are flagged as 
false detections, the latter are merged with the corresponding 
candidate in the union catalogue) or not associated (e.g., SZ 
sources clearly distinct from the MCXC clusters with no sig- 
nificant re-extracted signal at the cluster position and size). 

Finally, for MCXC clusters without redshift and luminosity 
information, the association was only based on distance, setting 
Dx < 5 arcmin, and the consolidated based on visual inspection 
of SZ, RASS and DSS images and other ancillary information. 
Two cases were found to be a mis-identification. The SZ candi- 
date was closer by chance to a faint XCS cluster, in the vicin- 
ity of the real counter-part (another MCXC cluster and an Abe 11 
cluster, respectively). 

4.3.2. Identification with optical clusters 

The Planck SZ candidates are associated with known clus- 
ters from optical catalogues (Abell, Zwicky, SDSS-based cat- 
alogues) on the basis of distance with a positional matching 
within a search radius set to five arcmin. The consolidation of 
the association was performed using the RASS information as 
described below, which allows us to mitigate the chance associ- 
ations with poor optical galaxy groups and clusters. 

SDSS-based catalogues - We have considered the four cat- 
alogues listed in Sect. 4.2. We define a quality criterion for the 
association, 2sdss, in terms of cluster richness as a proxy of 
the cluster mass (see for instance Johnston et al. 2007; Rozo 
et al. 2009). We set the quality criterion, 2sdss, to 3 for low 
reliability, to 2 for good reliability and to 1 for high relia- 
bility. The corresponding richness thresholds are 110 and 70 
for GsDSS = 1/2/3, respectively. The corresponding estimated 
masses (given the Msqq-Ri* relation) are M500 > 6.5 x 10^"^ M© 
and M500 > 3.5 X 10^"^ M©. However due to the large scatter 
and associated uncertainty in the mass estimate from the mass- 
richness relation, we consolidate the association of the Planck 
candidates with SDSS clusters by combining the 2sdss with the 
RASS signal at the Planck-csindidaiQ position (see Sect. 4.1.1). 
In practice, only associations with Qsdss = 1 or 2 and a signal- 
to-noise, measured at the Planck position in an aperture of five 
arcmin in the RASS survey, (S/N)rass ^ 1 are retained as firm 
identifications. We stress that our choice of richness thresholds 
is relatively conservative on average. Indeed, our 2sdss = 1 and 
2 matched candidates are found with high (S/N)rass values as 
shown in Fig. 17, with mean (S/N)rass = 7.1 and 6.6 and me- 
dian (S/N)rass = 5.9 and 5.4 for Qsdss = 1 and 2 matches, 
respectively. 
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Fig. 17: Normalized distribution of the signal-to-noise in RASS 
survey at the position of Planck SZ detections with SDSS 
richness-based quality 2sdss = 1 (solid line) and 2sdss = 2 
(dashed line). 



4.3.3. Identification with SZ clusters 

The association with known SZ clusters was performed within a 
five-arcmin radius. A visual inspection of the ancillary data and 
an a posteriori check of the RASS signal at the position of the 
Planck candidates associated with clusters from SZ catalogues 
is performed. It confirms that the values of (S/N)rass, when the 
coverage is significant, are high with an average value of 5.4. 

4.3.4. Identifications from NED and SIMBAD 

The information provided from querying NED and SIMBAD 
databases is mainly redundant with cross-checks with cluster 
catalogues. However, it lets us avoid missing a few associations. 
We therefore performed a systematic query in SIMBAD and 
NED with an adopted search radius set to five arcmin. Similarly 
to the association with clusters in optical catalogues, the posi- 
tional association is consolidated using the results of the search 
in RASS data. Furthermore, the Planck-csLndidaics solely match- 
ing NED or SIMBAD entries were inspected and the identifi- 
cation was confirmed or discarded using the information from 
WISE counter-parts and the DSS images. 



Abell and Zwicky catalogues - The Planck candidates are 
associated with Abell and Zwicky clusters on the basis of a 
positional matching within five arcmin. In the present case, 
we do not make use of any richness information in order 
to consolidate the association. We rather use here solely the 
RASS signal, (S/N)rass, at the SZ-candidate position. Planck- 
candidates associated with Abell or Zwicky clusters and with 
(S/N)rass ^ 1 are retained as firmly identified. For associations 
with (S/N)rass < 1 , we decided on a firm identification only af- 
ter checking the status of the counter-part in the WISE data and 
performing a visual inspection of the SZ signal and of the images 
from ancillary data, including DSS images. 



5. Follow-up programme for confirmation of 
Planck candidates 

We have undertaken, since Spring 2010, an extensive follow- 
up programme in order to perform a cluster-by-cluster confir- 
mation of the Planck cluster candidates and obtain a measure- 
ment of their redshifts. A total of 276 Planck candidates, se- 
lected down to S/N = 4 from intermediate versions of the Planck 
SZ catalogue, were observed in pursuit of their redshift mea- 
surement. We have constructed our strategy for the selection 
of the Planck targets primarily on the successful results of the 
series of follow-up observations in X-rays based on Director's 
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Table 4: List of the main observing facilities used for the confirmation of the Planck newly-discovered clusters, and for the mea- 
surement of their redshifts. 



Site 


Telescope 


Aperture (m) 


Instrument 


Filters 


Redshift 




XMM-Newton 




EPIC/MOS & PN 




Fe K 


La Palma 


NOT 


2.56 


ALFOSC 




spec 


La Palma 


INT 


2.5 


WFC 


griz 


phot 


La Palma 


GTC 


10.4 


OSIRIS 




spec 


La Palma 


TNG 


3.5 


DOLORES 




spec 


La Palma 


WHT 


4.2 


ACAM 


griz 


phot 


La Silla 


NTT 


3.7 


EF0SC2 




spec 


La Silla 


MPG/ESO 2.2m 


2.2 


WFI 


VRI 


phot 


Mullard Radio Ast. Obs. 


AMI 


3.7 and 13 


SA&LA 


13.5 to 18 GHz 




Tenerife 


IAC80 


0.82 


CAMELOT 


griz 


phot 


TUBITAK Nat. Obs. 


RTT 


1.5 


TFOSC 


gri 


spec,phot 



Discretionary Time on the XMM-Newton observatory (Planck 
Collaboration IX 201 1 ; Planck Collaboration Int. 1 2012; Planck 
Collaboration Int. IV 2013). Snapshot observations, sufficient to 
detect extended X-ray emission associated with Planck clusters 
and to estimate redshifts from the Fe line for the brightest clus- 
ters, were conducted sampling the SZ detections down to S/N 
= 4. These observations allowed us to better understand the SZ 
signal measured by Planck and hence to refine the criteria to se- 
lect targets, especially for further optical follow-up. 

We have engaged numerous campaigns on optical facilities, 
which now constitute our main means of confirmation of Planck 
SZ detections. Planck candidates with low-quality DSS images 
or without SDSS information, or low (S/N)rass, were primar- 
ily sent for deeper multi-band imaging observations. They were 
followed-up to the depth needed for the confirmation, i.e., find- 
ing an optical counter-part, and for the determination of a pho- 
tometric redshift. Candidates with galaxy concentrations in DSS 
or with counter-parts in SDSS, and/or with high (S/N)rass, were 
preferentially sent for spectroscopic confirmation. The priority 
being to confirm the clusters and to secure the largest number 
of robust redshifts, no systematic spectroscopic confirmation of 
photometric redshifts was performed for low-redshift clusters 
fephot < 0-4). For higher-redshift clusters, spectroscopic confir- 
mation of the photometric redshifts is more crucial. As a result, 
we have made use of telescopes of different sizes, from 1-m to 
10-m class telescopes, optimizing the selection of targets sent 
to the different observatories (Table 4 gives the list of the main 
telescopes). Eight- and ten-metre class telescopes, e.g., GTC, 
GEMINI and VLT, were used to spectroscopically confirm red- 
shifts above 0.5 for already confirmed clusters. 

Our efforts to confirm the Planck cluster candidates, mea- 
sure redshifts, and characterize cluster physical properties relies 
on ongoing follow-up of a large number of cluster candidates 
in the optical (ENO, RTT150, WFI), in the infrared (Spitzer^^) 
and at SZ wavelengths (Arcminute Microkelvin Imager, AMI). 
The output of the confirmation and redshift measurements from 
the observing campaigns is summarized in Sect. 6.2. Companion 
publications, in preparation, will detail the observing campaigns 
and their results. 

5. 1. XMM-Newton observatory 

The X-ray validation follow-up programme of 500 ks observa- 
tions undertaken in XMM-Newton DDT is detailed in Planck 



Under Spitzer programs 80162 and 90233. 



Collaboration IX (2011), Planck Collaboration Int. I (2012), 
and Planck Collaboration Int. IV (2013). It consisted of ob- 
serving 51 Planck targets and led to the confirmation of 43 
Planck cluster candidates, two triple systems and four double 
systems. There were eight false candidates. This follow-up pro- 
gramme has constituted the backbone of the Planck cluster con- 
firmation and most importantly has allowed us to better un- 
derstand the SZ signal measured by Planck and thus to bet- 
ter master the criteria for confirmation (or pre-confirmation) of 
the Planck cluster candidates. By providing us with the phys- 
ical properties and redshift estimates of the confirmed clus- 
ters, it has furthermore given us a first view on the phys- 
ical characteristics of the newly discovered Planck clusters. 
Snapshot observations (around 10 ks) of the Planck candidates 
took place between May 2010 and October 201 1. All the results 
from the four observing campaigns were published in Planck 
Collaboration IX (201 1), Planck Collaboration Int. I (2012), and 
Planck Collaboration Int. IV (2013). Calibrated event lists were 
produced with vll.O of XMM-Newton-S AS, and used to derive 
redshifts and global physical parameters for the confirmed clus- 
ters (Planck Collaboration IX 2011). The redshifts were esti- 
mated by fitting an absorbed redshifted thermal plasma model 
to the spectrum extracted within a circular region corresponding 
to the maximum X-ray detection significance. Most of the red- 
shifts were confirmed using optical observations. Additional ob- 
servations at VLT were conducted to confirm spectroscopically 
the highest redshifts. 

5.2. Optical observation in tine Northern hemispliere 
5.2.1. ENO telescopes 

In total 64 cluster candidates from Planck were observed at 
European Northern Observatory (ENO^^) telescopes, both for 
imaging (at IAC80, INT and WHT) and spectroscopy (at NOT, 
GTC, INT and TNG), between June 2010 and January 2013. 
The aims of these observations were the confirmation, photo- 



Observations are conducted under programme 090A-0925. 

ENO: http : / / www .iac.es/eno. php?lang=en. 

The observations were obtained as part of proposals for the Spanish 
CAT time (semesters 2010A, 2010B, 2011 A, 201 IB, 2012A and 
2012B), and an International Time Programme (ITP), accepted by the 
International Scientific Committee of the Roque de los Muchachos 
(ORM, La Palma) and Teide (OT, Tenerife) observatories (reference 
ITP12_2). 
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metric redshift measurement, and spectroscopic confirmation of 
redshifts above z = 0.3. 

INT, WHT and IAC80 - The optical imaging observations 
were taken either with the Wide-Field Camera (WFC) on the 
2.5-m Isaac Newton Telescope (INT), the auxiliary-port cam- 
era (ACAM) at the 4.2-m WilHam-Herschel Telescope (WHT), 
or with CAMELOT, the optical camera at the 0.82-m telescope 
(IAC80). The targets were observed in the Sloan gri filters. For 
the majority of fields, either Sloan z or Gunn Z images are also 
available. Images were reduced using the publicly-available soft- 
ware Iraf and Sextractor (Berlin & Amouts 1996). The 
data reduction included all standard steps, i.e., bias and flat field 
corrections, astrometric and photometric calibrations. The pho- 
tometric calibration is based either on standard star observations 
or, if available, on data from the SDSS. Finally, all magnitudes 
were corrected for interstellar extinction, based on the dust maps 
by Schlegel et al. (1998). We obtained photometric redshifts us- 
ing the BPZ code (Bemtez 2000), using a prior based on SDSS 
data, and fitting a set of galaxy templates. The BPZ code pro- 
vides the Bayesian posterior probability distribution function for 
the redshift of each object, which is later used in the process of 
cluster identification. The identification of the galaxy overden- 
sity located near the Planck positions and the estimate of the 
photometric redshifts of the associated clusters were performed 
using a modified version of the cluster-algorithm described in 
Sect. 4.1.2. 

GTC and TNG - Spectroscopic observations were performed 
using the 10-m Gran Telescopio Canarias (GTC) telescope and 
the 3.6-m Telescopio Nazionale Galileo (TNG) telescope. The 
OSIRIS spectrograph at GTC was used in long-slit mode to ob- 
serve a total of eight targets with two slit positions per candidate. 
We used the R500R grism and a binning 2x1, which provides a 
resolution R = 300 with a slit width 1 arcsec, and a wavelength 

coverage 4800-10000 A. We retrieved three exposures of 1200 s 
each. The final spectra present a S/N of about 20 in galaxies with 
/ = 20 mag. We used the DOLORES multi-object spectrograph 
(MOS) at TNG to observe 9 candidates. The masks were de- 
signed to contain more than 30 slitlets, 1.5 arcsec width, placed 
within an area about 6 arcmin x8 arcmin in order to cover the 
target field. We used the LR-B grism, which provides a disper- 
sion of 2.7 A/pixel, and a wavelength coverage between 4000 

and 8000 A. We carried out three acquisitions of 1 800 s each 
and obtained spectra with S/N ^ 15 in galaxies with / = 20 mag 
using a total integration time of 5400 s. 

Nordic Optical Telescope (NOT) - Spectroscopic redshift 
measurements were obtained using the Andalucia Faint Object 
Spectrograph and Camera (ALFOSC) at the NOT. Most targets 
were observed in MOS mode, targeting typically ten to fifteen 
galaxies per ALFOSC field (covering 6.4 x 6.4 arcmin^, with an 
image scale of 0.188 arcsec/pixel). One or two unfiltered 300s 
pre-imaging exposures were obtained per candidate cluster, in 
addition to a single 300s exposure in each of the SDSS g-and 
/ bands. The de-biased and flat field calibrated pre-imaging data 
were used to select spectroscopy targets. The final mask design 



The observing runs took place on June 28 - July 3, 2011, January 
20-25, 2012, July 16-21, 2012 and January 9-14, 2013. 

The MOS masks were cut at the Niels Bohr Institute, Copenhagen 
University. 



was carved out using custom software, generating slits of fixed 
width 1.5 arcsec and of length typically 15 arcsec. Grism No. 
5 of ALFOSC was used, covering a wavelength range 5000 - 

10250 A with a resolution of about R = 400 and dispersion 

3.1 A/pixel. Redwards of 7200 A strong fringing is present in 
the ALFOSC CCD. It was effectively suppressed using dither 
pattern alternating the placement of the spectroscopy targets be- 
tween these sets of slits. 

In addition to the MOS observations, spectroscopic observa- 
tions in single- slit mode were conducted for some Planck candi- 
dates. For these observations, a long slit covering the entire 6.4 
arcmin length of the ALFOSC field and a width of 1.3 arcsec 
was employed, with the same grism and wavelength coverage as 
for the MOS observations. The field angle was rotated to place 
the long slit over multiple targets, to include the apparent BCG 
as well as two to three other bright cluster galaxies within the 
ALFOSC field. 

5.2.2. RTT150 

A total of 88 Planck cluster candidates were followed up with 
the Russian Turkish Telescope (RTT150^^) from July 2011 
to December 2012 within the Russian quota of observational 
time. In total, about 50 dark nights, provided by Kazan Federal 
University and Space Research Institute (IKI, Moscow), were 
used for these observations. Direct images and spectroscopic 
redshift measurements were obtained using TUBITAK Faint 
Object Spectrograph and Camera (TFOSC^^), similar in layout 
to ALFOSC at NOT (see above) and to other instruments of this 
series. 

The TFOSC CCD detector cover a 13.3 x 13.3 arcmin^ area 
with 0.39 arcsec per pixel image scale. Direct images of clus- 
ter candidates were obtained in Sloan gri filters, in series of 
600s exposures with small 10-30 arcsec) shifts of the tele- 
scope pointing direction between the exposures. All standard 
CCD calibrations were applied using Iraf software, individ- 
ual images in each filter were then aligned and combined. The 
total of 1 800 s exposure time in each filter was typically obtained 
for each field, longer exposures were used for more distant clus- 
ter candidates. Deep multi-filter observations were obtained for 
all candidates, except those unambiguously detected in SDSS. 
With these data, galaxy clusters can be efficiently identified at 
redshifts up to z ~ 1 . 

Galaxy clusters were identified as enhancements of surface 
number density of galaxies with similar colours. Cluster red se- 
quences were then identified in the colour-magnitude diagram 
of galaxies near the optical centre of the identified cluster. The 
detected red sequence was used to identify the BCG and cluster 
member galaxies. Using the measured red-sequence colour pho- 
tometric redshift estimates were obtained, which were initially 
calibrated using the data on optical photometry for galaxy clus- 
ters from the 400SD X-ray galaxy cluster survey (Burenin et al. 
2007). 

For spectroscopy we used the long-slit mode of the instru- 
ment with grism No. 15, which covers the 3900-9100 A wave- 
length range with 12 A resolution when a slit of 1.8 arc- 
sec width is used. Galaxy redshifts were measured through the 
cross-correlation of obtained spectra with a template spectrum of 
an elliptical galaxy. Spectroscopic redshifts were typically ob- 



http : //he a .iki.rssi.ru/rttl50/en/index .php. 
2^ http: //hea. iki . rssi . ru/rttl50/en/ 
index . php?page=t f osc. 
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tained for the spectra of a few member galaxies, including the 
BCG, selected from their red sequence in the imaging observa- 
tions. These data allow us to efficiently measure spectroscopic 
redshifts for clusters up to z ~ 0.4. For the highest-redshift clus- 
ters, complementary spectroscopic observations were performed 
with the BTA 6-m telescope of SAO RAS using SCORPIO focal 
reducer and spectrometer (Afanasiev & Moiseev 2005). 

5.3. Optical observation in the Southern hemisphere 

5.3.1. MPG/ESO 2.2-m Telescope 

Optical imaging of 94 Planck cluster candidates in the 
Southern hemisphere was performed under MPG programmes 
at the MPG/ESO 2.2-m telescope using the Wide-Field Imager 
(WFI).^^ The WFI detector is a mosaic of 8 2k x 4k CCDs, cov- 
ering a total area of 33 arcmin x34 arcmin on the sky, with an 
image scale of 0.238 arcsec/pixel. Each field was observed in 
the V-, R-, and /-bands with a default exposure time of 1800s 
(with five dithered sub-exposures) per passband. The basic data 
calibration, including de-biasing and flat-field frame calibration, 
followed standard techniques. The individual exposures were re- 
registered and WCS calibrated using the USNO-B 1 catalogue as 
an astrometric reference before being stacked into a combined 
frame for each filter, covering the entire WFI field. Photometric 
redshifts of the observed clusters were then determined from an 
algorithm that searches for a spatial galaxy overdensity located 
near the position of the SZ cluster candidate that also corre- 
sponds to an overdensity mV - R versus R - I colour-colour 
space. The median colour of galaxies located in this overden- 
sity was then compared to predicted colours of early-type galax- 
ies at different redshifts by convolving a redshifted elliptical 
galaxy spectral energy distribution template with the combined 
filter+telescope+detector response function. 

5.3.2. New Technology Telescope (NTT) 

Observations^^ were conducted at the 3.5-m NTT at the ESO 
observatory at La Silla to measure spectroscopic redshifts of 33 
Planck clusters with the EF0SC2 instrument in the MOS mode. 
A clear BCG was identified in the clusters in pre-imaging data, 
and besides the BCG a redshift was measured for at least one 
other member of the cluster. In the following a brief outline of 
the observations and the data reduction are given (see Chon & 
Bohringer 2012, for details). 

Each field of the Planck target candidates was optically im- 
aged in Gunn r band for target selection and mask making. The 
imaging resolution is 0.12 x 0.12 arcsec^, and the field of view 
is 4.1 X 4.1 arcmin^ for both imaging and spectroscopic obser- 
vations. When necessary, the field was rotated to optimize target 
selection. We used the grism that covers the wavelength range 

between 4085 A and 7520 A, with 1.68 A per pixel at resolu- 
tion 13.65 A per arcsec. We typically applied 10 to 15 slitlets per 
field with a fixed width of 1 .5 arcsec for the MOS and of 2.0 arc- 
sec for the long-slit observations. Including at least three bright 



Based on observations under MPG programmes 086.A-9001, 
087.A- 9003, 088.A-9003, 089.A-9010, and 090.A-9010. The obser- 
vations were conducted during the periods of November 27 - December 

3, 2010, March 8-19, May 21 - June 3, and November 30 - December 

4, 2011, December 30, 2011 - January 7, 2012, June 10-18, 2012, and 
January 6-13 2013. 

The observations were performed during three spectroscopic ob- 
serving campaigns, 087.A-0740, 088.A-0268 and 089.A-0452. 



objects, preferably stars, to orient the field, the slitlets were allo- 
cated to the candidate member galaxies. The exposure times for 
the clusters range from 3600s to 10800s. 

The data were reduced with the standard reduction pipeline 
of I raf . The redshifts from the emission lines were determined 
separately after correlation with the passive galaxy templates. 
We use the rvsao package, which applies the cross-correlation 
technique to the input templates of galaxy spectra to measure the 
object redshift. The REFLEX templates were used for this analy- 
sis, which include 17 galaxy and stellar templates. We confirmed 
a spectroscopic cluster detection if at least three galaxies have 
their R- value greater than 5, and lie within ± 3000km/s of the 
mean velocity of the cluster members. We then took the median 
of those galaxy redshifts as the cluster redshift. For the long- 
slit observation, the cluster was confirmed with the redshift of 
the BCG and another galaxy at similar redshift within the afore- 
mentioned criteria. 



5.4. Observations in the SZ domain with AMI 

An ensemble of 60 Planck blind SZ candidates, spanning a range 
of S/N between 4 and 9 and meeting the Arcminute Microkelvin 
Imager (AMI) observability criteria, was observed with AMI. 
The goal of this programme was to confirm Planck cluster candi- 
dates through higher-resolution SZ measurements with AMI and 
to refine the position of confirmed clusters in order to optimize 
the subsequent optical follow-up observations aiming at redshift 
measurement. AMI comprises of two arrays: the Small Array 
(SA); and the Large Array (LA). Further details of the instru- 
ment are given in AMI Consortium et al. (2008). Observations 
carried out with the SA provide information that is well coupled 
to the angular scales of the SZ effect in clusters, whereas snap- 
shot observations obtained with the LA provide information on 
the discrete radio- source environment. The latter allowed us to 
detect the presence of nearby, bright radio sources, helping in 
further selecting the targets for observation with the SA. Details 
of the AMI data reduction pipeline and mapping are described 
in Planck and AMI Collaborations (2013). 



6. Results of the validation and follow-up 

The external validation allows us to identify Planck SZ detec- 
tions with previously-known clusters and to assemble crucial in- 
formation on the identified clusters such as their redshifts. The 
validation steps corresponding to the association with known 
clusters were performed following a chosen hierarchy: X-ray 
clusters from the updated MCXC meta-catalogue; then optical 
clusters from Abell and Zwicky catalogues; then optical clus- 
ters from the SDSS-based catalogues; followed by SZ clusters 
from SPT and ACT samples; and finally clusters from NED and 
SIMBAD queries. The first identifiers of the Planck SZ detec- 
tions given in Table B.l reflect the validation hierarchy. 

In the following, we present the results of the external val- 
idation process and of the follow-up campaigns for confirma- 
tion of Planck candidates and measurement of their redshifts 
(see Table 5 and Fig. 18). We also present the confirmation from 
SDSS galaxy catalogues and from X-ray archival data. We fur- 
ther discuss the unconfirmed candidate new clusters detected by 
Planck, which we classify into three categories of different reli- 
ability. 
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Table 5: Summary of the external validation and confirmation from follow-up observations. In each category, the number represents 
the total number of Planck candidates identified exclusively with X-ray, optical, or SZ clusters. The category X-ray clusters covers 
identifications from the updated MCXC meta-catalogue. The category optical clusters covers identification from the Abell, Zwicky, 
and SDSS-based catalogues only. The category SZ clusters covers identification from SPT and ACT catalogues only. Confirmations 
from follow-up do not cover the observations performed by the Planck collaboration to measure the missing redshifts of known 
clusters. Confirmation from archival data covers X-ray data from Chandra, XMM-Newton, and ROSAT PSPC pointed observations 
only. 
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Fig. 18: Distribution of the Planck clusters and candidates in the 
different categories defined in the external validation process. 
The validation follows the order association with MCXC clus- 
ters, then Abell and Zwicky clusters, then SDSS clusters, then 
SZ clusters, and finally clusters from NED/SIMBAD. 



6. 1. Planck clusters associated with known clusters 

A total of 683 out of 1227 SZ detections in the Planck cata- 
logue, i.e., 55.7%, are associated with previously-known clus- 
ters from X-ray, optical, or SZ catalogues, or with clusters found 
in the NED or SIMBAD databases. We give the number of clus- 



ters identified in each category and we discuss notable cases of 
known clusters that are not included in the Planck SZ catalogue. 



6.1 .1 . Identification with known X-ray clusters 

A total of 472 Planck SZ-candidates are identified with known 
X-ray clusters from the MCXC meta-catalogue, which repre- 
sents 38.5% of the Planck SZ detections and 69.1% of the iden- 
tifications with previously-known clusters. These identifications 
of course account for many Abell clusters in the RASS -based 
catalogues of X-ray clusters. 

Using the cluster properties reported in the MCXC and the 
Planck noise maps at the cluster positions, we computed the ex- 
pected SZ signal and the expected S/N ratio for a measurement 
with Planck. We have compared the number of detected clusters 
in the Planck catalogue with S/N > 4.5 to the number MCXC 
clusters at an expected significance of 4.5. Only 68 clusters ex- 
pected to be detected at S/N > 4.5 are not included in the Planck 
catalogue, including 16 with predicted S/N between 4 and 4.5. 
Of the 52 clusters with expected S/N > 4.5, only 41 are outside 
the masked regions and could thus be in the PSZ catalogue. Our 
computation of the expected SZ signal and S/N were based on 
scaling relations for X-ray- selected clusters, not accounting for 
the dispersion in the relations. We therefore focus on the non- 
detected MCXC clusters that significantly depart from the ex- 
pected S/N value, namely by more than 5 cr. A total of 13 clus- 
ters are in this category. The two objects RXCJ225 1.7-3206 and 
RXCJO 117.8-5455 show emission in high-resolution Chandra 
imaging that is point-like rather than extended and are likely not 
clusters of galaxies (Mantz et al. 2010; Magliocchetti & Briiggen 
2007). Of the other eleven missing MCXC clusters, some present 
AGN contamination. This is the case for RXC Jl 326.2+ 1230 
(Magliocchetti & Briiggen 2007), RXJ1532.9+3021 (Hlavacek- 
Larrondo et al. 2012), RXCJ1958.2-3011, RXCJ225 1.7-3206, 
and RXCJOl 17.8-5455 (Magliocchetti & Briiggen 2007), Abell 
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689 (Giles et al. 2012), ZwC12089 (Rawle et al. 2012), PKS 
0943-76 (Abdo et al. 2010), and Abell 2318 (Crawford et al. 
1999). In these cases, the presence of the AGN affects the 
X-ray luminosity measure leading to an overprediction of the 
SZ signal. Some exhibit significant radio contamination, e.g., 
RXCJ1253.6-3931 (Plagge et al. 2010) and RXCJ1958.2-3011 
(Magliocchetti & Briiggen 2007), which hampers the SZ de- 
tection. Cool-core clusters for which the X-ray luminosity is 
boosted due to the central density peak have an over-estimated 
expected SZ signal. This is the case for RXCJ0425. 8-0833 
(Hudson et al. 2010), ZwC12701 (Rawle et al. 2012), Abell 
1361 (Rafferty et al. 2008), and RBS 0540 (Eckert et al. 2011; 
Belsole et al. 2005). Other "missing" clusters are CIZA clusters: 
RXC J0643.4+4214, RXC J1925.3+3705, RXC J2042. 1+2426 
and RXC J0640.1-1253, REFLEX cluster RXCJ2149.9-1859, 
APMCC 699, Abell 3995, Abell 2064 and RBS 171. 

In addition to the clusters discussed above which are not in- 
cluded in the catalogue due to contamination by AGN or pres- 
ence of cool-cores etc., we note that some notable nearby ex- 
tended clusters are also not included in the Planck SZ catalogue. 
Indeed, the detection methods used to detect the SZ effect are 
not optimized for the detection of sources with scale radius ^500 
in excess of 30 arcmin. Of the 25 clusters in this category (with 
z < 0.03) in the MCXC meta-catalogue, six are included in the 
Planck catalogue. The remaining 19 fall into the masked areas 
(seven out of 19, among which Perseus and Abell 1060 lie in the 
PS mask (Fig. 19, first two panels), and Ophiuchus and 3C 129.1 
lie in the Galactic mask (Fig. 19, second two panels)) and/or 
have a S/N ratio below the PSZ catalogue threshold S/N = 4.5. 
This is the case of Virgo cluster (Fig. 19, lowest panel), which is 
detected in the Planck survey but with a signal-to-noise ratio at 
its position of about 3.9. Virgo's extension on the sky (^500 = 168 
arcmin) further hampers its blind detection. 

We show in Fig. 19 the reconstructed SZ signal from the 
MILCA algorithm (Hurler et al. 2010) for five of the "miss- 
ing" extended clusters. These clusters, despite not being part 
of the Planck catalogue of SZ sources, are well detected in the 
Planck survey. They all are included in the thermal SZ map con- 
structed from the Planck channel maps and presented in Planck 
Collaboration XXI (2013). 

6.1.2. Identification with known optical clusters 

A total of 182 Planck SZ detections are identified exclusively 
with optical clusters from Abell and Zwicky catalogues, and 
from the SDSS-based published catalogues, i.e., 26.6% of the 
known clusters in the Planck catalogue. 

The Planck SZ candidates at S/N > 4.5 have 111 exclu- 
sive associations with Abell or Zwicky clusters, i.e., with clus- 
ters not in any of the catalogues compiled in the MCXC meta- 
catalogue. In addition to these associations, 72 Planck detections 
are solely identified with clusters from the SDSS-based cata- 
logues. These are either rich and massive systems {Rl^ greater 
than 110, 2sDSS = 1 clusters) or moderately low-richness sys- 
tems (2sDSS = 2 clusters, exhibiting hot gas as indicated by 
their S/N value in the RASS survey). However, not all the rich 
GsDss = 1 clusters in SDSS-based catalogues are found in the 
Planck catalogue. A total of 213 2sdss = 1 clusters from all 
four SDSS-based catalogues (201 outside the Planck union PS 
and Galactic mask) are not included in the Planck catalogue. 

We explore why these rich clusters are not detected blindly 
by the SZ-finder algorithms. We first compare the richness-based 
masses against the X-ray luminosity-based masses of 26 of these 
"missing" clusters found in the MCXC meta-catalogue. We find 



a median ratio of 2.6 ±1.2 for the richness-to-X-ray based 
masses, indicating that the richness-based masses seem to be 
systematically overestimated. Unlike the X-ray clusters, we thus 
cannot compute a reliable estimate of the expected S/N value 
for SZ detection of these optical clusters. We therefore directly 
search for the SZ signal at the positions of the 201 "missing" 
SDSS-clusters and found that all of them have S/N values below 
the Planck threshold, with a mean signal-to-noise of 1.6, except 
for three clusters. Two of these three "missing" SDSS-clusters 
have their S/N value from the extraction at the cluster position 
slightly higher than 4.5. The increase in S/N value is due to the 
difference in estimated background noise when centring the ex- 
traction at the cluster position as opposed to the blind detection. 
The third missing rich cluster is affected by contamination from 
CMB anisotropy, which results in a bad estimate of its size and 
consequently of its SZ signal. 

6.1.3. Identification with known SZ clusters 

The majority of the SZ clusters, from SPT or ACT, used in the 
validation process are low-mass systems (M^^^^^^ around 2.3 x 
10^^ Mq). Planck is particularly sensitive to massive rich clusters 
and thus only a total of 56 of these clusters match Planck SZ 
detections, out of which 16 candidates are exclusively associated 
with SZ clusters^^ from ACT or SPT. Nine more ACT and SPT 
clusters are associated with Planck SZ detections between S/N = 
4 and 4.5. We have searched for the SZ signal in the Planck data 
at the position of the remaining non-observed ACT/SPT clusters 
by extracting the SZ signal at their positions. We found that all 
had signal-to-noise values lower than 4. 

We have also checked the redundancy of SZ detections 
within Planck by comparing the ESZ sample, constructed from 
10 months of survey with a cut at Galactic latitudes of ± 14 de- 
grees, with the present Planck catalogue. Of the 189 high sig- 
nificance ((S/N)esz > 6 ESZ detections, 184 ESZ confirmed 
clusters are included the present Planck catalogue within a dis- 
tance of five arcmin from their ESZ position. The mean sepa- 
ration between the ESZ and present positions is of order 1.35 
arcmin, within Planck'^ positional accuracy. Their S/N values 
were increased by a factor 1.17 on average with respect to their 
(S/N)esz, (Fig. 20) and only four out of six of the ESZ clus- 
ters have new S/N values significantly lower than ESZ signal- 
to-noise threshold (S/N)esz = 6. They are displayed as stars in 
Fig. 20. Four ESZ clusters are not included the present Planck 
catalogue, they fall in, or nearby, the PS mask used for the pre- 
processing of the channel maps prior to running the detection 
algorithms. Such a mask was not utilized for the construction of 
ESZ sample. We choose not to a posteriori include these four 
"missing" ESZ clusters in the present Planck SZ catalogue. 

6.1 .4. Identification with clusters from NED or SIMBAD 

As expected only a small number of clusters are identified from 
querying the databases, supplying identifiers for thirteen SZ 
Planck detections. This is because the information in NED and 
SIMBAD is redundant with that in the X-ray, optical, or SZ cat- 
alogues used for the external validation. The thirteen clusters 
found solely from querying the databases are found in the RASS 
survey but not in dedicated cluster catalogues, and thus not in- 
cluded in the MCXC; they are found in serendipitous Chandra 
surveys, or they are part of miscellaneous cluster catalogues. 

Six Planck clusters were confirmed from XMM-Newton or NTT ob- 
servations and are also published in Reichardt et al. (2013). 
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Fig. 19: Five selected nearby and extended clusters not included in the PSZ catalogue. All exhibit an extended SZ signal detected in 
the Planck survey. From top to bottom: Perseus cluster and Abell 1060 (in the point-source mask); Ophiuchus cluster and 3C 129.1 
(in the Galactic mask); and Virgo cluster (below the S/N threshold of the catalogue). Left panels: reconstructed thermal SZ maps 
from the MILCA algorithm (Hurler et al. 2010). The dashed circles represent the apertures of ^500 from the MCXC catalogue. 
Each SZ-map covers an area of 4^500 x 4^500. Right panels: composite images of the optical (DSS), X-ray (ROSAT) and SZ signal 
(Planck). The size of the composite images is adapted to optimize the display (Perseus: 2x2 square degrees; Abell 1060: 1 X2J 
square degree; Ophiuchus: 1x1 square degree; 3C 129.1: 0.77 x 0.77 square degrees; Virgo: 3.84 x 3.84 square degrees). The solid 
circle in the left comer shows, both in the SZ Planck map and in the composite image, a 10 arcmin field of view; except for Virgo 
for which it shows a 30 arcmin field of view. 
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Fig. 20: Ratio of the signal-to-noise in the present Planck cata- 
logue, (S/N)psz, to that in the ESZ sample (Planck Collaboration 
VIII 2011), (S/N)esz, for 184 confirmed ESZ clusters included 
in the Planck catalogue. Four clusters with signal-to-noise in 
the PSZ catalogue significantly smaller than the ESZ threshold 
((S/N)esz = 6) are shown as stars. 



6.2. Newly-discovered Planck clusters and candidates 

Among the 544 Planck SZ sources, we distinguish two cate- 
gories: (1) confirmed clusters, i.e., those that have been con- 
firmed by the follow-up programmes of the Planck collabora- 
tion^"^ or using the SDSS galaxy catalogues. We also add eight 
confirmations from X-ray archival data; (2) Candidate clusters 
with different levels of reliability, namely, CLASS 1 cluster candi- 
dates, that fulfil high-quality criteria for the SZ detection and for 
the associations and/or counterparts in ancillary data, CLASS 2 
candidate clusters, i.e., those that fulfil, on average, good-quality 
criteria, and CLASS 3, low-reliability cluster candidates. 



Confirmation from Planck collaboration follow-up pro- 
grammes At S/N > 4.5, a total of 233 Planck SZ detections 
were followed up in X-rays, optical, and SZ at the different fa- 
cilities listed previously, with some observations targeted to the 
measurement of spectroscopic redshifts for already known clus- 
ters. In total 157 Planck SZ detections with S/N > 4.5 were con- 
firmed as new clusters. Some of the Planck-confLrmcd clusters 
were also reported in recent cluster catalogues in the optical, e.g.. 
Wen et al. (2012) or in the SZ e.g., Reichardt et al. (2013). 
The analysis of the observations of Planck sources by AMI 
yielded ten sources with strong Bayesian evidences that have 
clearly visible decrements and were considered as confirmed, in- 
cluding the confirmation of three associations with optical clus- 
ters. 

For the candidates confirmed by XMM-Newton and by opti- 



A handful of new Planck clusters from the ESZ sample were con- 
firmed independently from the Planck collaboration by SPT (Story et al. 
2011), AMI (AMI Consortium et al. 2011), Bolocam (Sayers et al. 
2012a) and CARMA (Muchovej et al. 2012). 



cal telescopes, redshifts from Fe lines and from photometric or 
spectroscopic data are available. The validation of Planck clus- 
ter candidates with XMM-Newton has shown its particular effi- 
ciency in confirming SZ candidates due both to the high sensi- 
tivity of XMM-Newton, allowing Planck clusters to be detected 
up to the highest redshifts (Planck Collaboration XXVI 2011), 
and the tight relation between X-ray and SZ properties. The de- 
tection of extended XMM-Newton emission and a comparison 
between the X-ray and SZ flux permits an unambiguous confir- 
mation of the candidates. By contrast, confirmation in the optical 
may be hampered by the Planck positional accuracy and by the 
scatter between the optical observables and the SZ signal, which 
increase the chance of false associations. The XMM-Newton 
follow-up programme yielded 51 bona fide newly-discovered 
clusters, including four double systems and two triple systems. 
There were eight false candidates. Thirty-two of the 51 individ- 
ual clusters have high-quality redshift measurements from the 
Fe line. The relation between the X-ray and SZ properties was 
used to further constrain the redshift of the other clusters ; most 
of these redshifts were confirmed clusters using optical obser- 
vations. Out of a total of 37 single clusters confirmed by XMM- 
Newton, 34 are reported in the Planck catalogue of SZ sources 
at S/N > 4.5. Additionally four double systems are included in 
the present PSZ catalogue and were also confirmed by XMM- 
Newton. 

The follow-up observations conducted with optical telescopes 
lead to the confirmation and to the measurement of spectroscopic 
or photometric redshifts (companion publications, in prepara- 
tion, will present the detailed analysis and results from these 
follow-up). In the Northern hemisphere, 26 spectroscopic red- 
shifts for Planck clusters detected at S/N > 4.5 and observed 
at the RTT150 are reported, to date, in the PSZ catalogue. 
A dozen additional spectroscopic redshifts were measured for 
known clusters. Confirmation of 21 Planck SZ clusters detected 
above 4.5 were obtained with the ENO facilities (at INT, GTC 
and NOT), and robust redshift measurements were obtained 
for 19 of them, including 13 spectroscopic redshifts. In the 
Southern hemisphere, WFI observations provided photometric 
redshifts for 54 clusters included in the Planck catalogue at S/N 
> 4.5, while 19 spectroscopic redshifts obtained with the NTT- 
EF0CS2 instrument are reported in the Planck catalogue. 



Confirmation from SDSS galaxy catalogues The firm 
confirmation of the candidates was done through the follow-up 
observations for confirmation and measurement of their redshift 
as detailed above. However in the case of the Planck candidates 
falling in the SDSS footprint we also used the SDSS galaxy cat- 
alogues to search, as presented in Sect. 4.1.2, for galaxy over- 
densities associated with Planck SZ detections. This provides us 
with an estimate of the photometric redshifts, and in some cases 
we could retrieve spectroscopic redshifts for the BCG as well. 

In this process, the major uncertainty in the associations 
of Planck SZ detections with galaxy overdensities is due 
to chance associations with low-richness systems or associ- 
ations with diffuse concentrations of galaxies in the SDSS 
data. The XMM-Newton confirmation programmes (see Planck 
Collaboration Int. IV (2013) for discussion) showed that Planck 
candidates with SDSS counterparts were confirmed including 
PLCK G193.3-46.1 at z ^ 0.6. However, the X-ray analysis 
of the Planck detections with SDSS counterparts illustrated the 
difficulty in distinguishing between associations of Planck SZ 
signals with massive clusters and with pre-virialized structures. 
In particular, in the case of extended filamentary structures or 
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dynamically perturbed sources, an offset between the BCG po- 
sition and the concentration barycentre is noted. 

We considered the Planck SZ candidates with counterparts in 
the SDSS data taking into account diagnostics such as the rich- 
ness/mass estimates as well as the offsets between the SZ, the 
BCG and the barycentre positions. We further used the outputs 
of the search in WISE and in RASS data, and the associated im- 
ages, in order to assess the significance of the galaxy overdensity 
in SDSS at the position of the Planck candidates. For the Planck 
SZ detections where both ancillary data and SDSS barycen- 
tre/BCG positions agreed, we set that they are confirmed. We 
found a total of 13 such associations for which we report the 
photometric or the spectroscopic redshifts. It is worth noting 
that firm confirmation of these associations is needed and needs 
to be performed using either optical spectroscopic observations 
or X-ray observations of the Planck SZ detections. In the cases 
where the offsets between barycentre and BCG position output 
by the search in SDSS data were too large, and/or when other 
ancillary information was unable to discriminate between reli- 
able or chance associations, we have chosen to keep the status 
of candidate for the Planck SZ detection. These cases some- 
times also coincide with association of Planck detections with 
clusters from the SDSS cluster catalogues, with a quality flag 
GsDSS = 0, or with confusion in the association, i.e., with posi- 
tions not in agreement between counterpart and published SDSS 
clusters. We provide a note for all these cases in order to indicate 
that an overdensity in SDSS data was found. 

Candidate new clusters The remaining 366 Planck SZ 
sources, not identified with previously known cluster nor con- 
firmed by follow-up observation or ancillary data, are distributed 
over the whole sky (Fig. 21) and are yet to be firmly confirmed 
by multi- wavelength follow-up observations. They are charac- 
terized by an ensemble of quality flags defined in Sects. 4.1.1, 
4.1.2, and 4.1.3 based on the systematic searches for counter- 
parts in the public surveys during the external validation process. 
We further define an empirical Planck-'mitmdX quality flag Q^^. 
It assesses the reliability of the SZ detection itself from three in- 
dependent visual inspections of the nine Planck frequency maps, 
of frequency maps cleaned from Galactic emission and CMB, 
and of reconstructed j-maps or j-maps produced from compo- 
nent separation methods (e.g., Hurler et al. 2010; Remazeilles 
et al. 201 1). Moreover, we visualize the SZ spectra from the SZ- 
finder algorithms and from aperture photometry measurements 
at the candidate positions. Finally we correlate, at the position 
of the Planck SZ candidates and within an area of 10 arcmin 
radius, the j-map to the 857 GHz channel map, as a tracer of 
the dust emission, and to the Planck mono-frequency CO map at 
217 GHz (Planck Collaboration XIII 2013). The qualitative flag 
Q^^ combines all this information into three values 1 to 3 from 
highest to lowest reliability with the following criteria: 

- Q^^ = 1, i.e., high reliability: (i) Clear compact SZ source in 
the SZ maps; (ii) significant measurements of the SZ decre- 
ment below 217 GHz and good or reasonable detection at 
353 GHz; (iii) no correlation with dust nor CO emission and 
no rise of the 545 and 857 GHz fluxes on the thermal SZ 
spectrum. 

- = 2, i.e., good reliability: (i) visible SZ detection in 
the SZ map or significant detection of the SZ signal below 
217 GHz; (ii) contamination causing rise of the 545 GHz and 
possibly 857 GHz flux on the SZ spectrum without a strong 
correlation with dust and CO signals. 



- Q^^ = 3, i.e., low reliability: (i) weak SZ signal in the y- 
maps and/or noisy SZ maps; (ii) weak or no SZ signal in the 
cleaned frequency maps (iii) strong correlation (>80%) with 
dust and CO emission contamination with rising fluxes on 
the SZ spectrum at high frequencies, 353 GHz and above. 

We combine the qualitative SZ quality flag with the infor- 
mation from the search in the all-sky surveys, RASS and WISE, 
for counterparts of Planck candidates in order to assess the over- 
all reliability of the cluster candidates. We thus distinguish three 
classes of candidates: 

• CLASS 1 candidates. Highly -re liable candidates or pre- 
confirmed clusters: these are the Planck SZ detections that 
have a high probability of being associated with bona fide 
clusters and need to fulfil high-quality criteria for SZ, RASS, 
and WISE detections. We retain in this category Planck SZ 
detections with high or good SZ quality flags (Q^^ = 1 or 
2) and with a RASS-BSC source (not coinciding with stars) 
or with (S/N)rass ^ 2, i.e., SZ detections with quality flag 
Grass = 1- The CLASS 1 candidates furthermore have to ful- 
fil a condition of high or good probability (>80%) of being 
associated with an overdensity of galaxies in the WISE sur- 
vey. 

We find 54 CLASS 1 Planck candidates ranging from S/N of 
4.5 to 6.3, with a median signal-to-noise ratio of 4.8. The 
majority of them are detected by two methods and 25.9% of 
them are detected only by one method. They are distributed 
as 26 and 28 Q^^ = 1 and 2 candidates, respectively. These 
candidates show significant X-ray emissions with a median 
(S/N)rass - 3.7 and a mean of 4.2. 

• CLASS2 candidates. Reliable cluster candidates: they repre- 
sent 170 Planck SZ detections that show good or high quality 
criteria either in SZ or in RASS or in WISE without fulfilling 
all of them at once. Amongst them 61 have = 1 and 109 
have Q^^ = 2. 

• CLASS3 candidates. Low-reliability cluster candidates: 
these Planck SZ detections are the poor-quality, = 3, 
detections. They can also be associated with good quality, 

= 2, detections for which there are no good indications 
of the presence of an X-ray counterpart ((S/N)rass < 0-5 
and high probability of false association with ESC sources 
>2.5%) or a counterpart in the WISE survey (probability of 
association <70%). 

This class of candidates contains 142 Planck SZ detections 
with 27 and 115 SZ detection of quality Q^^ = 2 and 3, 
respectively. 

It is worth noting that this definition of the CLASS 3 Planck can- 
didates is dominated by the assessment of the SZ quality com- 
plemented by information from ancillary data. In doing so we 
assemble in this category of candidates the SZ detections that are 
either false or very low quality due to contamination. Moreover, 
according to the statistical characterization from simulations, 
about 200 false detections are expected. The number of false 
detections could be smaller since the simulations do not repro- 
duce the entire validation procedure, in particular omitting the 
cleaning from obvious false detections. Figure 22 suggests that 
the CLASS 3 candidates are likely to be dominated by false detec- 
tions. Therefore, we would like to warn against dismissing entire 
CLASS 3 of the catalogue as populated with false detections as 
some CLASS 3 candidates may be real clusters. For this reason, 
we choose not to remove these detections from the PSZ cata- 
logue but rather flag them as low-reliability candidates. Careful 
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Fig. 21: Distribution of the Planck SZ candidates across the sky. Blue symbols represent the CLASS 1 candidate clusters and red the 
CLASS 2 candidates. The open symbols stand for the CLASS 3 low-reliability SZ sources. 



follow-up programmes are needed in order to separate real clus- 
ters of galaxies from false detections among the CLASS 2 and 
CLASS3 objects. 

In order to illustrate our classification defined in terms of re- 
liability, we stack the signal in patches of 2.51 degrees across, 
centred at the position of the Planck clusters and candidates 
in the nine channel maps of Planck, removing a mean signal 
estimated in the outer regions where no SZ signal is expected 
(see Fig. 22 with the rows arranged from 30 GHz, upper row, to 
857 GHz, lower row). The stacked and smoothed images are dis- 
played for the Planck SZ detections identified with known clus- 
ters, CLASS 1, CLASS2 and CLASS3 candidates. Fig. 22 from left 
to right column. We clearly see the significant detection of both 
the decrement and increment of the 683 Planck clusters and of 
the Planck candidates of CLASS 1 and CLASS2. For the Planck 
SZ detections associated with bona fide clusters the increment is 
clearly seen at 353 and 545 GHz and is detected at 857 GHz. The 
smaller sample of the CLASS 1 highly reliable candidates shows, 
in addition to the decrement at low frequency, a good detection 
of the increment at 353 GHz. The significance of the increment 
at 545 GHz is marginal and no signal is seen at 857 GHz. The 
case of the CLASS 2 candidates (good reliability) shows that we 
now have lower-quality SZ detections (62% of the CLASS 2 can- 
didates have a good but not high SZ quality flag). This is illus- 
trated by the fact that an excess emission is detected at 217 GHz, 
most likely due to contamination by IR sources, and both at 545 
and 857 GHz where emission from dust is dominating. As for 
the stacked signal of the CLASS 3 sample of low-reliability can- 
didates, it does not show any significant SZ detection across fre- 
quencies, as compared to the sample of Planck detections iden- 
tified with known clusters (Fig. 22, right column). This confirms 
on statistical grounds the definition of the sample dominated by 



definition by the low-quality SZ, Q^^ = 3, detections repre- 
senting 84% of the detections in this class. Not surprisingly, the 
stacked signal of the CLASS 3 candidates shows a large amount of 
contamination across all Planck frequencies. The low-frequency 
signal is dominated by radio contamination, and/or CO emission 
at 100 GHz, while the high-frequency signal is contaminated by 
emission from dust or extragalactic point sources. A more quan- 
titative analysis is presented in Sect. 7.1. 

6.3. Summary of the external validation and redshift 
assembly 

The Planck catalogue of SZ sources comprises a total of 861 
identified or confirmed clusters with only nine percent of them 
being detected by one SZ-finder algorithm. We summarize in 
Table 5 and Fig. 18 the results of the cluster identification. 
Figure 23 illustrates the status of the Planck SZ detections. In 
particular, 70.2% of the Planck SZ detections with S/N>4.5 have 
so far been associated with clusters. The fraction increases to 
about 73% at S/N = 6. 

We have assembled, at the date of submission, a total of 
813 redshifts for the 861 identified or confirmed Planck clus- 
ters, which we provide together with the published Planck cat- 
alogue. Their distribution is shown in Fig. 24. In the process of 
the redshift assembly that is summarized below, especially for 
the already known clusters, we have favoured homogeneity for 
the sources of redshift rather than a cluster-by-cluster assembly 
of the most accurate z measure. A large fraction of the redshifts, 
456 of them, shown as the dashed green histogram in Fig. 24 
correspond to the spectroscopic redshifts quoted in the updated 
MCXC meta-catalogue (Piffaretti et al. 2011). They are associ- 
ated with the Planck clusters identified with known X-ray clus- 
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Fig. 22: Stacked signal in the nine Planck frequencies (30 to 857 GHz from upper to lower row). From left to right are displayed 
the Planck SZ detections identified with known clusters, the CLASS 1 high-reliability Planck SZ candidates, the CLASS 2 good- 
reliability Planck SZ candidates, and finally the CLASS 3 low-reliability SZ sources. The three lowest- frequency- channel images 
were convolved with a 10 arcmin FWHM Gaussian kernel, whereas the remaining six highest- frequency-channel images were 
smoothed with a 7 arcmin FWHM Gaussian kernel. 
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Fig. 23: Status of the Planck SZ sources. Left-hand-axis plots 
show the distribution of all Planck sources (in red). The blue 
line represents the known or new confirmed clusters and, among 
these, the clusters with a reported redshift measurement in black. 
Right-hand-axis cumulative distributions show, as a function of 
signal-to-noise, the fraction of known or new confirmed clusters 
in blue and those with a redshift in black. 



ters and they are denoted Planck-MCXC. For the Planck-MCXC 
clusters without reported redshifts from the MCXC, we have 
complemented the information with the available redshifts from 
NED and S 1MB AD. We have further quoted when available, 
mainly for the MACS clusters, the estimated photometric red- 
shifts from SDSS cluster catalogue of Wen et al. (2012). At the 
end only two Planck detections identified with MCXC clusters 
remain without redshifts. The redshift distribution of the Planck 
clusters identified with MCXC clusters mostly reflects that of the 
REFLEX/NORAS catalogues at low and moderate redshifts and 
the MACS clusters at higher redshifts. 

For the Planck detections exclusively identified with Abell or 
Zwicky clusters, we choose to report the redshifts published 
in the NED and SIMBAD data bases rather than those quoted 
in the native catalogues. As for the Planck detections identi- 
fied with clusters from the SDSS-based catalogues, we choose 
to favour homogeneity by reporting whenever possible the Wen 
et al. (2012) redshifts. Furthermore, we favour when available 
spectroscopic redshifts over photometric ones. The Planck de- 
tections exclusively associated with ACT or SPT clusters have 
pubhshed redshifts (Sifon et al. 2012; Hasselfield et al. 2013; 
Reichardt et al. 2013). We select in priority the spectroscopic 
ones when available. If not, we quote the photometric redshifts. 
Finally, the follow-up observations for confirmation of Planck 
detections started in 2010 and are still ongoing. As mentioned 
earlier our priority was to assemble the largest possible num- 
ber of confirmations and redshifts. Therefore, we did not sys- 
tematically confirm the photometric redshift estimates spectro- 
scopically. We report the obtained redshifts when available. In 
some cases, the new Planck clusters were confirmed from imag- 
ing or pre-imaging observations and the analysis is still ongo- 
ing. The spectroscopic redshifts will be updated when available. 
Spectroscopic redshifts for some known clusters will also be up- 
dated. A dozen Planck clusters were confirmed by a search in the 
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Fig. 24: Distribution of redshifts for the Planck SZ clusters 
(black line). The Planck clusters associated with MCXC clusters 
are shown in dashed green and the brand new Planck clusters are 
in the filled red histogram. 

SDSS galaxy catalogues. For these clusters, only a photometric 
redshift estimated by the cluster-finder algorithm of Fromenteau 
et al. (2013) is available and is reported. 

We show in Fig. 24 the distribution of redshifts of the Planck 
clusters. The mean redshift of the sample is 0.25 and its median 
is 0.22. One third of the Planck clusters with measured redshifts 
lie above z = 0.3. The new Planck clusters probe higher red- 
shifts and represent 40% of the z > 0.3 clusters. Their mean 
redshift is 0.38 and the median is z = 0.35. At even higher red- 
shifts, z > 0.5, the Planck catalogue contains 65 clusters includ- 
ing Planck SZ clusters identified with WHL12's clusters (Wen 
et al. 2012), or with clusters from ACT and SPT, or with X- 
ray clusters. The Planck detections in this range of redshifts, 29 
Planck new clusters, almost double the number of high redshift 
clusters. 

The Planck SZ catalogue has been followed up by the Planck 
collaboration using different facilities and only a small fraction 
of the Planck candidates were observed to date. A systematic 
follow-up effort for the confirmation of the remaining cluster 
candidates will likely reveal clusters at redshifts above 0.3. As a 
matter of fact, very few new clusters were found below z = 0.2 
(see Fig. 24). Such an observational programme is challenging 
and will most likely be undertaken by the Planck collaboration 
and by the community. It will increase further the value of the 
Planck SZ catalogue as the first all- sky SZ- selected catalogue. 

7. Physical properties of Planck SZ clusters 

The first goal of the external validation process based on the an- 
cillary multi- wavelength data is to assess the status of the Planck 
SZ detections in terms of known clusters, brand new clusters 
or cluster candidates. The wealth of information assembled and 
used during this process also allows us to explore the proper- 
ties of the Planck SZ clusters and candidates. We present in the 
following some of these properties, namely the contamination 
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levels of the Planck SZ detections, a refined measurement of the 
Compton Y parameter for the Planck clusters identified with X- 
ray clusters from the MCXC, an SZ-mass estimate based on a 
new proxy for all the Planck clusters with measured redshifts, 
and an estimate of the X-ray flux from the RASS data for the 
Planck SZ detections not included in the X-ray catalogues. This 
additional information associated with the Planck clusters and 
candidates derived from the validation process is summarized in 
the form of an ensemble of outputs given in Table C.l . 

We further present an updated and extended study of the 
SZ versus X-ray scaling relation, confirming at higher precision 
the strong agreement between the SZ and X-ray measurements 
(within 7?5oo) of the intra-cluster gas properties found by Planck 
Collaboration XI (2011). 

7. 1. Point-source contamination 

Galactic and extragalactic sources, emitting in the radio or in- 
frared domain, are known to lie in galaxy clusters and hence are 
a possible source of contamination for the SZ measurement (e.g., 
Rubino-Martm & Sunyaev 2003; Aghanim et al. 2005; Lin et al. 
2009). We address the possible contamination of the SZ flux by 
bright radio sources that may affect the measured signal in the 
direction of some of the Planck SZ detections. In order to do so, 
we searched for known radio sources in the vicinity of the Planck 
cluster candidates. In particular, we use the NVSS 1.4 GHz sur- 
vey (Condon et al. 1998) and SUMMS 0.85 GHz survey (Bock 
et al. 1999) to identify bright radio sources within seven arcmin 
of the Planck cluster or candidate position. We assumed a spec- 
tral index a = -0.5 for these sources to extrapolate their flux 
to the Planck frequencies. Most bright sources in NVSS and 
SUMSS have steeper spectral indexes (-0.6 or -0.7), so the 
value a = -0.5 provides us with an upper limit in most cases. 
After convolving the radio sources by Planck's beam, we esti- 
mate the maximum amplitude in units of jjK within five arcmin 
of the Planck position. We report only those cases where this 
amplitude is above 5 jiK in the 143 GHz channel and could thus 
contaminate the SZ signal. Below this value, the emission from 
radio sources can be considered negligible. 

We find that a total of 274 Planck clusters and candidates, 
i.e., 22% of the SZ detections, are affected by such emission 
from bright radio sources. These clusters or candidates are iden- 
tified in the PSZ catalogue and a specific note is provided. We 
find that the fraction of contaminated Planck SZ clusters identi- 
fied with known X-ray, optical, or SZ clusters is also 22%. The 
Planck candidate-clusters of CLASS 1 and CLASS2 are less con- 
taminated by bright radio sources; only a fraction of 15% and 
17% for CLASS 1 and 2, respectively. This is due to the defini- 
tion of our quality criteria for SZ detection, which results in less 
contamination for the high and good reliability candidates. 

Another approach used to assess the contamination is based 
on the stacking analysis of the Planck clusters and candidates 
described in Sect. 6.2. This analysis is performed on the sam- 
ple of Planck clusters identified with known clusters and on 
the sample of low-reliability CLASS 3 Planck candidates. To do 
so we fit a GNFW pressure profile to the signal at 100 GHz 
and 143 GHz and we subtract the associated SZ signal from 
the stacked maps. The residual signal is then compared with a 
toy model for point sources (Fy = S ^^^iy/ 30 GHz)^''''^ for radio 
sources) and (Fy = 5'g^^(y/857 GHz)^^'^ for IR point sources). 
Note that the residual signal at high frequencies is a combina- 
tion of possible IR sources and IR emission from Galactic dust; 
the latter is not explicitly modelled in the present analysis. The 
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Fig. 25: Stacked spectrum for known clusters SZ fluxes across 
Planck frequency bands. Stacked fluxes are measured in an aper- 
ture equal to the FWHM of the 143 GHz channel (i.e., about 
7 arcmin) for the known clusters (black filled circles) and the 
low-reliability CLASS 3 candidates (red filled triangles). The as- 
sociated uncertainties correspond to the fluctuation of the back- 
ground outside the cluster region. The average signal is esti- 
mated in each channel before (upper panel) and after (lower 
panel) the removal of the SZ signal. The average signals ex- 
pected from IR and radio sources are shown as solid and dashed 
lines, respectively. Red and black lines are for CLASS 3 and bona 
fide clusters, respectively. No subtraction of an SZ signal is per- 
formed for the CLASS3 candidates. 



PS toy models are convolved by the beam at each frequency and 
the signal is measured at a fixed aperture set to the FWHM of the 
143 GHz channel. The average signal within this aperture is es- 
timated for each channel before (Fig. 25, upper panel) and after 
(Fig. 25, lower panel) removal of the SZ signal. The black filled 
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Fig. 26: Comparison of the different Y estimates for the Planck 
clusters identified with MCXC clusters. In green are the blind 
measured Y values and in black are the refined y5oo,PSX measured 
fixing the size and positions to the X-ray values. Both are plotted 
as a function of the new proxy Y^. 



circles are for Planck SZ sources associated with known clusters 
and the red filled triangles stand for CLASS 3 candidates. The av- 
erage signal from the PS models is shown in Fig. 25 as solid 
(IR sources) and dashed (radio sources) lines. Red and black are 
for CLASS 3 and bona fide clusters, respectively. The error bars 
correspond to the fluctuation of the background outside the clus- 
ter region. For the sample of CLASS 3 candidates no SZ-signal 
removal was applied, since no significant detection is seen at 
100 GHz or 143 GHz. 

We find that the residual signal (after SZ subtraction) in the 
sample of known Planck clusters is compatible with the emis- 
sion from radio sources at low frequencies with {S^^^^^a^^^) = 
(14.6 mJy, -1) for the known clusters. It is also compatible with 
IR emission at high frequencies with a spectral index aiR = 2.5, 
in agreement with the results of Planck Collaboration Int. VII 
(2013) and with = 0.117 Jy. For CLASS3, where no SZ sig- 
nal is subtracted, it is the full signal that is compatible with the IR 
emission at high frequencies, with (5'g^^,QfiR) = (43.9 Jy, 2.5), 
and with radio emission from point sources with (S^^^^arad) = 
(117.1 mJy, -0.8). 

7.2. Refined measurement of Y 

While the true Fsoo is expected to be a low-scatter mass proxy, 
this is not the case for the blind Y^qq. Without a cluster-size es- 
timate, Y500 cannot be accurately measured. Moreover, the blind 
SZ flux is biased high on average, because the size is over- 
estimated on average. This effect is amplified by the non-linear 
nature of the size-flux degeneracy, with a larger effect of size 
over-estimation than size under-estimation. This behaviour, first 
identified and discussed in Planck Collaboration VIII (201 1) and 
Planck and AMI Collaborations (2013), hampers the direct use 
of the blind SZ fluxes as a mass proxy. As shown in Planck 
Collaboration VIII (2011), this degeneracy calls for a refined 



measurement of the SZ signal. In this section, we present two 
ways of refining the Y measurement. Both are based on fixing 
the cluster size in two cases, by setting it equal to the X-ray esti- 
mated size or by using the redshift information when available. 
The outputs of the refined measurement are provided as addi- 
tional information complementary to the catalogue of Planck SZ 
detections (see Appendix C and Table C.l). 

7.2.1 . Y at fixed X-ray size and position 

As shown by Planck Collaboration VIII (201 1), the size-flux de- 
generacy can be broken by introducing a higher-quality estimate 
of the cluster size ^500. This prior is directly provided by X-ray 
observations using an X-ray mass proxy such as Ix or the lumi- 
nosity Lx. Resorting to estimates of the cluster size from opti- 
cal richness is also possible, but suffers from the large scatter in 
richness-mass relation, as discussed previously. 

A detailed investigation of the effects of fixing the clus- 
ter size was presented in Planck Collaboration XI (2011, 
Appendix A). Following this approach, and for the Planck de- 
tections identified with clusters from the MCXC meta-catalogue, 
we have adopted the R500 and z values reported in Piffaretti et al. 
(2011) as priors to re-extract at the X-ray position the SZ sig- 
nal denoted F5oo,psx assuming the Arnaud et al. (2010) pres- 
sure profile (see Table C.l). The comparison between the blind 
Y500 and refined F5oo,psx (Fig. 26) shows that both the scat- 
ter and the offset are significantly reduced by the refined SZ 
measure. The SZ re-extraction at X-ray position and fixing the 
size to the X-ray derived size provides an unbiased estimate of 
the SZ signal. However, as stressed in Planck Collaboration XI 
(201 1, Appendix A), the MCXC cluster size derivation involves 
the M5oo-^x,50o relation, which exhibits a non-negligible scat- 
ter. This leads to a remaining systematic discrepancy between 
the expected Y value from X-ray measurements and the actual 
SZ flux derived from the Planck data. The use of the Yx proxy 
does not suffer from such an effect, but high-quality X-ray data 
permitting the use of such a quantity are not available for a large 
number of clusters (see Sect. 7.5 for the presentation of a sample 
of Planck SZ clusters with high-quality X-ray data). 

7.2.2. Ffrom the Y(0)-M relation 

The size-flux degeneracy can further be broken, as proposed by 
Arnaud et al. (2013), using the M^oq-D^Ysoo relation itself that 
relates ^500 and Y500, when z is known. Then Ysoo is derived from 
the intersection of the Msoq-D^Y^oo relation and the size-flux 
degeneracy curve. A detailed description of the method and the 
comparison of results in terms of bias and scatter can be found 
in Arnaud et al. (2013). 

The derived Y50Q parameter is denoted Y^ (since it involves 
a measurement of the Compton Y signal for clusters with mea- 
sured redshift z). It is the SZ mass proxy Y^ that is equivalent 
to the X-ray mass proxy Yx- Y^ is computed for all the 813 
Planck clusters with measured redshifts. We use Malmquist- 
bias-corrected scaling relation between mass and Y given in 
Planck Coflaboration XX (2013) 
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with E^(z) = ^lm(l + z)^ + computed in the fiducial ACDM 
cosmology. 

In Fig. 26, the refined ^500 value, measured fixing the size 
and position to the X-ray values y5oo,PSX, is compared to the 
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Fig. 27: Distribution of masses for the Planck SZ clusters, known 
or new confirmed clusters (solid black line), compared to the 
distribution of masses from the RASS-based cluster catalogues 
(dashed blue line). The masses for the MCXC clusters are es- 
timated from the luminosity-mass relation. The masses for the 
Planck clusters are computed using the SZ-proxy. The filled red 
histogram shows the distribution of the newly-discovered Planck 
clusters. 



blind y as a function of the derived proxy. We see that the 
scatter and the offset are significantly reduced. 

Under the two hypotheses of cosmology and scaling rela- 
tion, Yz provides the best estimate of Ysoo for the Planck SZ 
clusters and conversely a homogeneously-defined estimate of 
an SZ-mass, X-ray calibrated, denoted mJ^^. For the ensem- 
ble of Planck clusters with measured redshifts, the largest such 
sample of SZ- selected clusters, we show in Fig. 27 the distri- 
bution (black solid line) of the masses obtained from the SZ- 
based mass proxy. The distribution of the SZ masses is compared 
with those of the RASS clusters (dashed blue line) computed 
from the X-ray luminosity Lx,5oo- The mean and median masses 
of the Planck clusters are 3.3 and 3.5 x 10^^ M©, respectively. 
The Planck SZ catalogues contains all the massive clusters of 
the RASS catalogues. Interestingly, the distribution of newly- 
discovered Planck clusters extends to higher masses with a me- 
dian mass of 5.7 x 10^^ M©. Besides providing a homogeneous 
estimate of the masses from an SZ proxy for the largest SZ se- 
lected sample of clusters, we show that Planck detections sig- 
nificantly extend the mass range in the high-mass region up to 
1.6 X 10^5 Mo. 

7.3. distribution and comparison with other surveys 

Based on the masses derived from the SZ-proxy, we illustrate 
for MMF3 the M-z distribution of Planck SZ clusters detected 
over 83.7% of the sky. We show in all panels of Fig. 28 the lim- 
iting mass Miim computed following Planck Collaboration XX 
(2013) for three values of the completeness: 20% (solid line); 
50% (dashed line); and 80% (dotted line). The upper left panel 
exhibits the Planck clusters, with redshifts, detected by MMF3 at 
S/N>4.5. The mass limit corresponds to the average limit com- 



puted from the noise over the 83.7% sky fraction used by the SZ- 
finder algorithm. The resulting Mum is not representative of the 
inhomogeneity of the noise across the sky (see Fig. 3). We there- 
fore show the limiting mass in three areas of the sky (Fig. 3): the 
deep-survey area (upper right panel); the medium-deep survey 
area (lower left panel); and the shallow- survey area (lower right 
panel). The lines indicate the limit at which clusters have C% 
chances to be detected (C being the completeness value). We 
clearly see that whereas the average at 20% completeness 
does not fully represent the SZ detections by MMF 3 , the limiting 
masses in different survey depths are more representative of the 
detection process. We further note that except at low redshifts, 
z < 0.3 - 0.4, the Planck cluster distribution exhibits a nearly 
redshift-independent mass limit with a cut that varies according 
to the survey depth. 

It is worth examining the distribution of the Planck SZ clus- 
ters in the M-z plane and comparing it to that of other cata- 
logues. For illustration, we compare to an X-ray selected sam- 
ple, namely REFLEX-I, on the one hand (Fig. 29, right panel 
green open circles) and to the large-area SZ-selected cluster cat- 
alogues by ACT (Hasselfield et al. 2013) and SPT (Reichardt 
et al. 2013), on the other hand (Fig. 29, red open symbols). In 
this comparison we report, for the ACT clusters (open squares), 
the so-called UP? (Universal Pressure Profile) masses given in 
Hasselfield et al. (2013). 

The range of redshifts covered by the Planck SZ sample, 
from z = 0.01 to about 1 with 67% of the clusters lying be- 
low z = 0.3, is quite complementary to the high redshift range 
explored by ACT and in by SPT. For the comparison of the mass 
distribution we take advantage of our newly-proposed SZ-mass 
estimate, derived from F^, which provides us with a homoge- 
neous definition of the masses over the whole range of Planck 
SZ clusters with measured redshifts. The Planck clusters pop- 
ulate the full redshift range and they quite nicely fill a unique 
space of massive, M > 5 x 10^^ M©, and high redshift z > 0.5 
clusters, as shown in Fig. 29. This contrasts with the SZ clusters 
detected in 720 square degrees of SPT observations and those 
of ACT observations, which are dominated, as shown in Fig. 29 
left panel, by lower-mass higher-redshift clusters (up to z ~ 1.3). 
The combination of Planck and SPT/ACT catalogues samples 
the M-z space in a complementary manner. Clearly the all- sky 
nature of the Planck makes the most massive clusters preferen- 
tially accessible to Planck whereas the highest redshift clusters, 
z> 1 , are accessible to SPT. 

Very few massive high-redshift clusters exist in the X-ray 
catalogues, as seen in Fig. 29 (right panel open blue squares). 
The all-sky NORAS/REFLEX catalogues (Bohringer et al. 
2000, 2004) are limited to z = 0.45, a result of the (1 -h z)"^ sur- 
face brightness dependence of the X-ray detection limit (Fig. 29, 
right panel solid green line). The smaller- area MACS sample, 
based on systematic follow-up of RO SAT bright sources (Ebeling 
et al. 2007), contains a dozen clusters at z > 0.5. The 400SD 
sample (Burenin et al. 2007), based on serendipitous detec- 
tions in 400 deg^ of ROSAT pointed observations, contains only 
two clusters with M > 5 x 10^"^ M© and z > 0.5. Finally, 
only a couple of clusters in the range M > 5 x 10^"^ M© are 
found in the XMM-Newton based serendipitous cluster samples 
(XCS, Mehrtens et al. (2012); XMM-LSS, Pacaud et al. (2007); 
XDCP, Fassbender et al. (2011)). By contrast to an X-ray se- 
lected cluster catalogue, the Planck detection-limit, illustrated 
for the medium-deep survey zone and shown in Fig. 29 (right 
panel solid black line), has a much shallower dependence on 
redshift and is quasi-redshift independent above z = 0.4. The 
difference in cluster selection starts at redshifts z > 0.2. As a re- 
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Fig. 28: Mass limit illustrated for SZ detections by MMF3 algorithm. Upper left: average mass limit computed from the average 
noise over the sky. Upper right: same for the deep survey zone corresponding to 2.7% sky coverage centred at the Ecliptic polar 
regions. Lower left: same for the medium-deep survey area covering 41.3% of the sky. Lower right: same for the shallow- survey 
area covering 56% of the sky. In each panel, only detections in the corresponding areas are plotted. The lines dotted, dashed and 
solid lines show the Planck mass limit at 80, 50 and 20% completeness, respectively. 



suit of the quasi-redshift independent mass-selection of SZ sur- 
veys, Planck probes deeper than the X-ray selection. This is also 
seen in the overall distribution of redshifts of the Planck clusters. 
Fig. 24. 

This leaves the Planck SZ catalogue as the deepest all- sky 
catalogue spanning the broadest cluster mass range from 0.1 to 
1.6 X 10^^ Mq, and particularly adapted to the detection of rare 
very massive clusters in the tail of the distribution in the range 
M>5x\0^^MQ^n&z>0.5. 

7.4. X-ray flux of the Planck clusters and candidates 

For all Planck SZ detections, we estimated the unabsorbed fluxes 
at Earth in the [0. 1-2.4] keV band (as in the MCXC) measured 



in an aperture of five arcmin. The aperture is centred on the 
Planck candidate position, except for candidates associated with 
a BSC source, for which we adopt the X-ray position, since the 
BSC source is very likely the counterpart (Planck Collaboration 
Int. IV 2013). The conversion between the RASS count rate in 
the hard band and flux is performed using an absorbed thermal 
emission model with the A^h value fixed to the 21 cm value. 
The conversion depends weakly on temperature and redshift 
and we assumed typical values of kT = 6keV and z = 0.5. 
Planck Collaboration Int. IV (2013) compared such flux esti- 
mates with precise XMM-Newton fluxes measured within R500, 
S 500, for candidates confirmed with the XMM-Newton follow-up 
programme. These clusters lie in the range 0.1 < z < 0.9 and 
the 0.3 X 10-^2 < 5 500 < 6 X 10"^^ erg s"^ cm'^ flux range. The 
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Fig. 29: Left panel: distribution in the M-z plane of the Planck clusters (filled circles) compared with the SPT clusters (open light red 
circles) from Reichardt et al. (2013) and ACT catalogue (open red squares) from Hasselfield et al. (2013). Right panel: distribution 
in the M-z plane of the Planck clusters (black symbols) as compared to the clusters from the REFLEX catalogue (green open circles) 
Bohringer et al. (2004). The black crosses indicate the Planck clusters in the REFLEX area. The open blue squares represent clusters 
from the MCXC catalogue with redshifts above z = 0.5. The green solid line shows the REFLEX detection limit whereas the black 
solid line shows the Planck mass limit for the medium-deep survey zone at 20% completeness. 




Fig. 30: Ratio between RASS flux, computed in an aperture of five arcmin in radius centred on the Planck position, and MCXC value 
for Planck candidates identified with MCXC clusters. The fluxes are computed in the [0.1-2.4] keV band at Earth and corrected for 
absorption. S 500 is the flux corresponding to the luminosity within 7?5oo published in the MCXC catalogue. Left panel: the ratio is 
plotted as a function of distance between the Planck and X-ray positions; Middle panel: same, as a function of cluster redshift, for 
distances smaller than five arcmin; Right panel: same as middle panel, for RASS flux within /?5oo derived from the aperture flux, 
using the MCXC iterative procedure based on the L500-M500 relation and the R EXCESS gas density profile (Piffaretti et al. 2011). 
The red line is the median ratio in distance or redshift bins with the grey area corresponding to ±1 cr standard deviation in each bin. 



RASS aperture fluxes were found to underestimate the "true" 
flux by about 30%. 

Figure 30 extends this comparison further to all the Planck 
SZ detections identified with MCXC clusters. Piffaretti et al. 
(2011) published homogenized L500 and 7?5oo values derived 
from the flux given in the original catalogues in various aper- 
tures, using an iterative procedure based on the R EXCESS 
L500-M500 relation and gas density profile shape. We simply 
computed S 500 from L500, taking into account the K-correction at 
the cluster redshift, but neglecting its variation with temperature. 



Although derived from ROSAT survey data as our present 
flux estimate, ^'500 values from the MCXC are expected to be 
more accurate due to: (i) optimum choice of the X-ray cen- 
tre; (ii) higher S/N detection; (iii) more sophisticated flux ex- 
traction adapted to data quality and source extent (e.g., growth 
curve analysis); and (iv) use of R^m rather than a fixed aperture. 
Not surprisingly, the ratio between the present flux estimate and 
the MCXC value decreases with increasing offset between the 
Planck position and X-ray position (Fig. 30, left panel). The ratio 
drops dramatically when the distance is larger than five arcmin, 
i.e., when the X-ray peak lies outside the integration aperture. 



33 



Planck Collaboration: Planck catalogue of Sunyaev-Zeldovich sources 




0.001 0.010 
[arcmin'] 



Fig. 31: X-ray unabsorbed flux versus SZ flux. For Planck SZ 
detections identified with MCXC clusters (open green circles), 
the X-ray flux is estimated from L500. For other Planck SZ detec- 
tions, the flux is derived from RASS count-rate in a five-arcmin 
aperture (see Sect. 4.1.1). Planck new clusters and CLASS 1 can- 
didates are shown as open red circles and squares, respectively. 
The two lines corresponds to the expected Lsoo-i^soo relation 
(Amaud et al. 2010) at z = 0.01 and z = 1, respectively. 



tified with MCXC clusters have low (S/N)rass values) and the 
relatively large scatter (±30% standard deviation). For z < 0.1 
clusters, and if the RASS detection is reasonably good a more 
precise procedure is recommended, such as an adapted growth 
curve analysis, on a case-by-case basis. 

7.5. Scaling relations between SZ and X-ray quantities 

A fundamental scaling relation is that between Fsoo and its X- 
ray analogue. Fx- Introduced by Kravtsov et al. (2006), Fx is the 
product of Mg,5oo, the gas mass within /?5oo, and Tx, the spectro- 
scopic temperature outside the core . From the fact that the gas 
density profile used to compute Mg^soo is derived from deprojec- 
tion of the X-ray surface brightness profile, and that the X-ray 
emission depends on the square of the density, the ratio of these 
two quantities is 

Cxsz^x QMrsooTx 

Q ^ V(^e)dr 
<^e)dr 

where the angle brackets denote volume-averaged quantities, 
and Q is the dumpiness factor at the scale of the radial bins used 
to derive the density profile. The numerical constant Cxsz = 
o-T/(mecVemp) = 1 .416 X lO'^^ Mpc^ (Mq keV)-^ The ratio 
thus depends only on the internal structure of the intra-cluster 
medium. 

The properties of the Fx-^^soo relation, in particular its vari- 
ation with mass and redshift and the dispersion about the mean 
relation, are important probes of the physics of cluster formation. 



Those are rare cases, 18 nearby clusters (z < 0.1 with a median 
value of z = 0.05), for which a physical offset likely contributes 
to the overall offset. When these cases are excluded, the median 
ratio is 0.85 and depends on redshift (Fig. 30, middle panel); it 
significantly decreases with decreasing redshift below z of 0. 1 . 
The median ratio is 0.65 and 0.92, with a standard deviation of 
0.10 and 0.15 dex, below and above z = 0.1, respectively. This 
is mostly due to the choice of a fixed aperture that becomes too 
small as compared to R500 at low z. If we apply the same iterative 
procedure used by Piffaretti et al. (2011) to estimate 5*500 from 
the aperture flux, the resulting value is consistent on average 
with the MCXC value at all redshifts (Fig. 30, right panel). The 
dispersion is slightly increased. The aperture unabsorbed fluxes 
are thus reliable estimates of the X-ray fluxes above z > 0.1 on 
average. 

Figure 3 1 shows the X-ray flux as function of Fsoo for Planck 
candidates identified with known clusters, for the confirmed new 
Planck clusters and for the CLASS 1 candidates. For Planck de- 
tections identified with MCXC clusters we plot the more precise 
published S 500 value. All three categories of sources behave in 
a similar manner in good agreement with the range of redshifts 
probed by the sample. In this respect CLASS 1 candidates do not 
exhibit any departure with respect to the known or confirmed 
clusters. We provide the X-ray fluxes for the Planck clusters 
and candidates that are not identified with MCXC clusters (see 
Appendix C and Table C.l). For the Planck cluster with MCXC 
identifier, we refer the reader to the RASS catalogue outputs or 
to the homogenized MCXC meta-catalogue. The main limitation 
of the aperture unabsorbed fluxes is the statistical precision on 
the RASS estimate (most of the Planck SZ detections not iden- 



7.5.1. Data set 

Here we extend the study of a sample of 62 clusters from the 
Planck-ESZ sample with good quality XMM-Newton archive 
data presented in Planck Collaboration XI (2011, hereafter 
PEPXI). This study found Z)^ Ysoo/Cxsz Yx = 0.95 ± 0.03, in 
a good agreement with REXCESS prediction, 0.924 ± 0.004, of 
Amaud et al. (2010). 

All 62 objects in the PEPXI sample are included in the 
present catalogue. We further add 40 clusters from the cata- 
logue, including nine additional objects from the XMM-Newton 
archival study of Planck-dctcctcd LoCuSS systems presented 
by Planck Collaboration Int. Ill (2013), and the 31 Planck- 
discovered clusters with good redshift estimates (Qz = 2) con- 
firmed with the XMM-Newton (Planck Collaboration IX 2011; 
Planck Collaboration Int. I 2012; Planck Collaboration Int. IV 
2013). The total sample thus consists of 102 clusters. 

For each object, Yx and the corresponding R500 value were 
estimated simultaneously by iteration about the M^qq-Yx rela- 
tion of Amaud et al. (2010), 



E^I\z)Msm = 10 



14.567 



2x lOi^MokeV 



0.561 



Mr: 



(7) 



In the present study, we focus on the physical relation between 
Fsoo and Fx- While these quantities must be estimated within the 
same radii, the exact value of R500 is irrelevant as the radial de- 
pendence of the Fsoo/Fx ratio is negligible. We thus propagated 
only the measurement uncertainties on the temperature and gas 



Here we use the temperature measured in the [0. 15-0.75] Rsoo aper- 
ture. 
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Fig. 32: Relation between the Comptonization parameters F500, and the normalized Fx parameter for a sub- sample of the present 
catalogue. Black points show clusters in the Planck-ESZ sample with XMM-Newton archival data presented by Planck Collaboration 
XI (2011) and additional LoCuSS clusters studied by Planck Collaboration Int. Ill (2013). Green points represent new Planck 
clusters confirmed with XMM-Newton (Planck Collaboration IX 2011; Planck Collaboration Int. I 2012; Planck Collaboration Int. 
IV 2013). The red line denotes the scaling relations of Planck Collaboration XI (201 1). Left panel: relation in units of arcmin^ where 
Fsoo is extracted using the Arnaud et al. (2010) pressure profile. The grey area corresponds to median Y^qq values in Yx bins with 
±cr standard deviation. Right panel: scaling relation between the intrinsic Compton parameter, D^Fsoo, and Fx for the sub-sample 
of S/N > 7 clusters used in the cosmological analysis. The data are corrected for Malmquist bias, and Ysoo is extracted using the 
Planck pressure profile (see text). The black line is the best-fit power-law relation. 



Table 6: The Ysoq-Y^ relation. Column(l): sample under consideration; Column(2): Malmquist bias correction; Column (3): Pressure profile shape 
used in F500 extraction; Column (4): number of clusters in the sample; Column (5-6) Slope and normalization of the best-fit relation expressed 
as YsoolY^ - A{YxIY^Y, using BCES orthogonal regression. The pivot is Fp = lO'^^Mpc^; Column (7-8) Intrinsic and raw scatter around the 
best-fit relation; Column (9): mean ratio in logarithm, A(log(2))=log(75oo/Fx); Column (10-11): corresponding intrinsic scatter and raw scatter. 
The scatters are error-weighted values. The best estimate is indicated in bold face. The last line gives the REXCESS prediction (Arnaud et al. 
2010). 
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mass profiles, fixing the aperture to R^qq. We ignored the statis- 
tical and systematic uncertainties on the Msoo-i^soo relation it- 
self. Similarly Fsoo was re-extracted at the X-ray position with 
size fixed to X-ray size. Its uncertainty corresponds to the statis- 
tical error on the SZ signal. The results are summarized Table 6, 
with the best estimate indicated in bold face. 



These must however be taken into account when using 7500 or Fx 
as a mass proxy, e.g., when calibrating the ^500-^500 relation from 
combining the M^oq-Y^ relation and the relation between F500 and 7x 
(or equivalently M500). This calibration is extensively addressed in the 
Planck Collaboration XX (2013). 



7.5.2. The best-fit Y^^-Yx relation 

The Isoo-i^c scaling relation for the full sample is shown in units 
of arcmin^ in Fig. 32. At high flux the points follow the PEP XI 
relation. The slope and normalization are determined at slightly 
higher precision, due to the better quality SZ data. The derived 
intrinsic scatter (Table 6) is significantly smaller, a consequence 
of the propagation of gas mass profile errors in the Fx error bud- 
get, which was neglected in our earlier study. 

The relation levels off at around Fx = 5 x 10"^ arcmin^, with 
a bin average deviation increasing with decreasing Yx (Fig. 32 
left panel). This is an indication of Malmquist bias, as noted by 
Planck Collaboration Int. I (2012). Full correction of this bias 
when fitting scaling relations involves drawing mock catalogues 
according to the cluster mass function, to which the sample se- 
lection criteria are then applied. The present sample is a small 
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subset of the full S/N > 4.5 Planck catalogue and thus such a 
procedure cannot be applied. To minimize bias effects we will 
only consider high S/N detections, S/N > 7. To correct for the 
residual bias, we adapted the approach proposed by Vikhlinin 
et al. (2009). Before fitting the Y^m-^x relation, each individual 
Y value was divided by the mean bias, /?, given by 

exp(-x^/2(T^) 

ln/^= , ^cr, (8) 

^[7^|2tr^c[x| ^Jlcr) 

where x = - logCy/Fmin), 5^min being the flux threshold corre- 
sponding to the signal-to-noise cut, (S/N)cut- At the location of 
the cluster, F/Fmin = (S/N)/(S/N)cut. Here cr is the log-normal 
dispersion at fixed Fx- We took into account both the intrin- 
sic dispersion (Tint, estimated iteratively, and the statistical dis- 
persion, given by o- = ^In [((S/N) -h 1)/(S/N)]^ -h [In lOo-intjl 
The correction decreases the effective Y^qq values at a given Fx, 
an effect that is larger for clusters closer to the S/N threshold; 
i.e., low-flux objects. The net effect on the scaling relation is 
small, giving a 0.7 cr decrease of the normalization and a slight 
steepening of the power-law slope (Table 6). 

The slope and normalization of the relation are robust to the 
inclusion of newly-discovered Planck clusters. The results de- 
rived from the extended sample of 78 clusters with S/N > 7 
agree with those obtained for the updated XMM-ESZ sample 
within 0.5 (T (Table 6). They are also in agreement with the sub- 
sample of 7 1 S/N > 7 clusters included in the cosmological sam- 
ple discussed by Planck Collaboration XX (2013). We measured 
a significant intrinsic scatter of (Tint = 0-07 ± 0.01 dex. There 
is one spectacular outlier with an Y500 / Yx ratio nearly twice as 
big as the mean. This is the Planck ESZ cluster identified with 
A2813 or RXC J0043.4-2037 in the REFLEX catalogue, located 
at z = 0.29. Its high ratio is very puzzling. It cannot result from 
an inaccurate redshift measurement, as this is based on spectro- 
scopic data for several cluster galaxies (Bohringer et al. 2004). 
There is no evidence of a peculiar dynamical state from the X- 
ray morphology, and there is no evidence of contamination in 
the SZ data. 

Part of the dispersion could be due to the use of an inap- 
propriate fixed pressure profile in the Y500 extraction. When in- 
cluding possible errors on Y500 due to dispersion around the 
mean Amaud et al. (2010) profile, the scatter is decreased to 
(Tint = 0.06, a decrease at the 1 cr level. To further assess the ef- 
fect of the choice of the pressure profile, we re-extracted the SZ 
signal using the Planck+XMM-Newton profile shape measured 
for ESZ clusters by Planck Collaboration Int. V (2013, here- 
after PIPV). Individual profiles are used for Planck ESZ clusters, 
and the mean profile is used for the other clusters. This should 
give the most reliable estimate of the Fsoo-Fx relation, since it is 
based directly on measured profile shapes. In this case, the slope 
and scatter remain unchanged but the normalization is slightly 
decreased (at the 0.5 (T level). This is a result of the more in- 
flated nature of the PIPV profile as compared to the Arnaud et al. 
(2010) R EXCESS profile. The relation derived using PIPV pres- 
sure profiles is plotted in the right-hand panel of Fig. 32 together 
with the corresponding data points. 

The relation does not exhibit significant evidence of variance 
of the Fsoo / Fx ratio with mass, the slope is consistent with unity, 
as expected from strong self- similarity of pressure profile shape. 
However, we found an intrinsic scatter about three times larger 
than the results of Kay et al. (2012). Partly this is due to the 
presence of outliers in our data set (as discussed above), or it 
may be due to projection effects in observed data sets (Kay et al. 



2012). The mean ratio is very well constrained with a precision 
of 2.5%,log(F5oo/Fx) = -0.027 ±0.010. This confirms at higher 
precision the strong agreement between the SZ and X-ray mea- 
surements (within 7?5oo) of the intra-cluster gas properties found 
by PEP XI and other studies (Andersson et al. 201 1 ; Sifon et al. 
2012; Marrone et al. 2012; Rozo et al. 2012). The ratio is per- 
fectly consistent with the X-ray prediction, suggesting that the 
dumpiness must be low. However there are still large systemat- 
ics that are discussed in Appendix D. 

8. Summary 

Planck'^ all-sky coverage and broad frequency range are de- 
signed to detect the SZ signal of galaxy clusters across the sky. 
We provide, from the first 15.5 months of observations, the 
largest ensemble of SZ-selected sources detected from an all- 
sky survey. The Planck catalogue of SZ sources contains 1227 
detections. This catalogue, statistically characterized in terms of 
completeness and statistical reliability, was validated using ex- 
ternal X-ray and optical/NIR data, alongside a multi-frequency 
follow-up programme for confirmation. A total of 861 SZ de- 
tections are confirmed associations with bona fide clusters, of 
which 178 are brand-new clusters. The remaining cluster candi- 
dates are divided into three classes according to their reliability, 
i.e., the quality of evidence that they are likely to be bona fide 
clusters. 

A total of 813 Planck clusters have measured redshifts rang- 
ing from z = 0.01 to order one, with one-third of the clusters 
lying above z = 0.3. The brand-new Planck clusters extend the 
redshift range above z = 0.3. For all the Planck clusters with 
measured redshift, a mass can be estimated from the Compton 
F measure. We provide a homogeneous mass estimate ranging 
from (0.1 to 1.6) x 10^^ M©. Except at low redshifts, the Planck 
cluster distribution exhibits a nearly redshift-independent mass 
limit and occupies a unique region in the M-z space of massive, 
M > 5 X 10^"^ Mo, and high-redshift (z > 0.5) clusters. Owing 
to its all-sky nature, Planck detects new clusters in a region of 
the mass-redshift plane that is sparsely populated by the RASS 
catalogues. It detects the rarest clusters, i.e., the most massive 
clusters at high redshift in the exponential tail of the cluster mass 
function that are the most useful clusters for cosmological stud- 
ies. With the presently confirmed Planck SZ detections, Planck 
doubles the number of massive clusters above redshift 0.5, as 
compared to other surveys. The Planck SZ catalogue is, and will 
be for years to come, the deepest all-sky SZ catalogue spanning 
the broadest cluster mass range. 

The Planck SZ catalogue should motivate multi- wavelength 
follow-up efforts. The confirmation of the cluster candidates 
will reveal clusters at higher redshifts than the present distribu- 
tion. Such follow-up efforts will further enhance the value of 
the Planck SZ catalogue as the first all- sky SZ selected cata- 
logue. It will serve as a reference for studies of cluster physics 
(e.g., galaxy properties versus intra-cluster gas physics, dynam- 
ical state, evolution, etc.). Using an extended sub-sample of 
the Planck SZ clusters with high-quality XMM-Newton data, 
the scaling relations between SZ and X-ray properties were re- 
assessed and updated. With better-quality data and thus higher 
precision, we show excellent agreement between SZ and X-ray 
measurements of the intra-cluster gas properties. We have thus 
derived a new up-to-date reference calibrated local relation be- 
tween F and Fx. 

The Planck SZ catalogue will also serve to define samples for 
cosmological studies. A first step in this direction is already 
taken in Planck Collaboration XX (2013), where an analysis of 
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the SZ cluster abundance to constrain the cosmological parame- 
ters is performed using a sub- sample selected from the PSZ cata- 
logue consisting of 189 clusters detected above a signal-to-noise 
ratio of 7 with measured redshifts. The value-added information 
derived from the validation of the Planck SZ detections, in par- 
ticular the SZ-based mass estimate, increases even further the 
value of the Planck SZ catalogue. 

The combination of the Planck all- sky SZ data with near 
future and planned observations of the large-scale structure by 
surveys such as PAN-STARRS, LOFAR, EucHd, LSST, and 
RSG/e-ROSITA will revolutionize our understanding of large- 
scale structure formation and evolution. 
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Appendix A: Selection of Frequency Channel 
Maps 

An assessment of which combination of Planck frequency chan- 
nels to use was performed using the MMF 1 implementation of the 



matched multi-filter described in Sect. 2.2.2. The HFI and LFI 
channel maps were preprocessed as described in Sect. 2.1, with 
the only difference being that the point- source mask contained, 
in addition, detections from the LFI channel maps with S/N 
> 10. Five different combinations of frequency channels were 
investigated, all Planck channels (30-857 GHz), all HFI channels 
plus the 70 GHz channel map from LFI (70-857 GHz), all HFI 
channels (100-857 GHz), the five lowest frequency HFI chan- 
nels (100-545 GHz) and the four lowest frequency HFI channels 
(100-353 GHz). For each combination of frequency channels a 
catalogue of SZ sources was extracted, resulting in five different 
catalogues; the only differences between them must be entirely 
due to the choice of channels in the combination. 

The first four of these catalogues are in good agreement in 
terms of the clusters detected, with all the differences amongst 
them being due to detections with S/N < 5. The (100-353 GHz) 
catalogue, however, contains significantly more detections, re- 
sulting in a poor agreement between it and the other catalogues 
that is not limited to low S/N detections. This is interpreted as 
being due to the lack of a dust-dominated channel in this com- 
bination, without which it is more difficult to constrain contam- 
ination due to dust emission. 

In order to assess any improvement in the S/N ratios of de- 
tected clusters with the inclusion of extra data, a robust sam- 
ple of reliable sources is required. To produce this, only clus- 
ters outside the 65% dust mask and with S/N > 8 were kept 
from each combination. The differences in the S/N of the same 
sources detected using different frequency channel combina- 
tions can then be examined. The ratio between the S/N values 
of the common detections in each combination to those of the 
(100-857) combination was then found; the mean of this ratio is 
shown in Table A.l. This approach clearly shows the (100-353) 
combination to be considerably noisier than the other combina- 
tions, which is consistent with the observations reported above. 
Neither the inclusion of the LFI frequency channels or just the 
70 GHz channel brings any significant improvement in the S/N 
of the clusters. Using the six HFI channel combination results 
in marginally better S/N than the (100-545 GHz) combination. 
The frequency channel combination chosen therefore is (100- 
857 GHz) since this gives the highest S/N with the smallest data- 
set. Reducing the S/N threshold from 8 to 6 and hence doubling 
the number of SZ sources used to evaluate the mean ratio does 
not change the conclusions of this analysis. 

Appendix B: Extract from the Planck catalogue 
of SZ sources 

We describe here the Planck catalogue of SZ sources delivered 
by the collaboration and available together with the individual 
lists from all three detections methods, the union mask used by 
these methods and the ensemble of notes on individual clusters 
at http : // www .sciops.esa.int/ index . php?page= 
Planck_Legacy_Archive&pro ject=planck . 

The union Planck SZ catalogue contains the coordinates 
and the signal-to-noise ratio of the SZ detections and a sum- 
mary of the validation information, including external identifi- 
cation of the cluster and redshifts if they are available. The ex- 
ternal identification quoted in the delivered product corresponds 
to the first identifier as defined in the external validation hier- 
archy, namely identification with MCXC clusters followed by 
Abell and Zwicky, followed by SDSS -based catalogues, fol- 
lowed by SZ catalogues, followed finally by searches in NED 
and SIMBAD. Due to the size-flux degeneracy discussed in 
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Table A. 1 : Mean of the ratio of the S/N of the common detections between the catalogues produced using different frequency 
channel combinations, excluding all detections inside the 35% dust mask. The improvement in the S/N of the detected clusters 
between the (100-353) and (100-857) combinations is clearly demonstrated, as is the lack of significant improvement in S/N when 
LFI data is included. The improvement between the (100-545) and (100-857) combinations is smaller and in the region of 1 to 2%. 



Combination 


(100-353)/(100-857) 


(100-545)/(100-857) 


(70-857)/(100-857) 


(30-857)/(100-857) 


S/N> 6 


0.86 


0.99 


1.00 


1.00 


S/N> 8 


0.84 


0.98 


1.00 


1.00 



Sect. 2.3, no reference flux quantity is outputted for the union 
catalogue. 



The individual catalogues from the three detection methods, 
MMFl, MMF3, and PwS, contain the coordinates and the S/N ra- 
tio of the detections, and information on the size and flux of 
the clusters. The size is given in terms of and the flux is 
given in terms of the total integrated Comptonization parame- 
ter, Y = y5R5QQ . The full information on the degeneracy between 
and Y is provided in the form of the two-dimensional marginal 
probability distribution for each cluster. 

The degeneracy information is provided in this form so that it 
can be combined with a model or external data to produce tighter 
constraints on the parameters. For example, combining it with 
an X-ray determination of the size can be done by taking a slice 
through the distribution at the appropriate • This is what is done 
in Sect. 7.2.1 and the refined measurement using X-ray informa- 
tion can be found in Table C.l . 



Table B.l presents an extract of the PSZ catalogue, in terms 
of the first rows of the online table and the following selected 
columns: 

NAME: name of cluster. 

R.A., DEC: right ascension (J2000) and declination (J2000). 
S/N: signal-to-noise ratio of the SZ detection. 
VALIDATION: status of the SZ detection from external 
validation: 20 = previously-known cluster; 10 = new confirmed 
Planck cluster; 1 = CLASS 1 candidate ; 2 = CLASS2 candidate; 
3 = CLASS 3 candidate. 

IDext* first external identifier of the known clusters, 
z: redshift of the cluster as reported from the external validation. 
COMMENTS: F = no comment; T = comment. Comments are 
readable in an external file. 

The complete version of the PSZ catalogue also contains the 
additional columns: 

INDEX: index of the detection, determined by the order of 
the clusters in the union catalogue and sorted into order of 
ascending Galactic longitude. 
GLON, GLAT: Galactic coordinates. 
POS_ERR: errors on the position. 

PIPELINE: pipeline from which information is taken; namely 
1 =MMF1; 2 = MMF3; 3 = PwS. 

PIPE_DET: pipeline making the detection, with the following 

order in bits: 1st = MMFl; 2nd = MMF3; 3rd = PwS. 

PCCS: flag for a match with sources from the PCCS catalogue. 

COSMO: flag for those clusters that are included in the sample 

used for the cosmological analysis of Planck Collaboration XX 

(2013). 



Appendix C: Outstanding outputs from the 
external validation 

Based on the ancillary data used for the validation of the Planck 
SZ catalogue, we provide value-added information to the Planck 
SZ detections. 

Namely, we provide, in addition to the first external identifier, 
possible other common identifiers, IDs. 

We report the redshift information associated with the Planck 
clusters (z) and specify its source, (scr). 

For clusters with measured redshifts, we compute the SZ-proxy 
Yz and the mass estimate (M^qo^ associated errors. For the 
clusters identified with MCXC clusters we provide the SZ signal, 
i^500,psx, re-extracted fixing the size to the X-ray size provided 
in the MCXC catalogue at the X-ray position. Note that the X- 
ray positions used in the present study are those quoted in the 
MCXC meta-catalogue. The positions reported in the ESZ sam- 
ple were taken from a sampled grid of coordinates with a pixel 
size of 1.71 arcmin. Due to this sampling, the reported MCXC 
positions in the ESZ sample exhibit an average offset of 70 arc- 
sec (less than a pixel, which varies depending on the position of 
the object on the sphere). 

For Planck SZ detections not associated with a previously- 
known X-ray cluster and with a signal-to-noise, (S/N)rass > 
1 (T, we provide the unabsorbed X-ray flux, Sx (and error), mea- 
sured in an aperture of 5 arcmin in the band [0.1-2.4] keV. We 
only provide an upper limit in the case of (S/N)rass < 1 cr, ex- 
cept for three SZ detections for which RASS exposure is very 
low and (S/N)rass < -5 cr. The aperture is centred on the Planck 
position, except for candidates associated with a BSC source for 
which we adopt the X-ray position. These clusters are flagged. 

Appendix D: Systematic effects on the X-ray 
versus SZ scaling relation 

Both X-ray and SZ measurements are likely affected by system- 
atic effects linked to e.g., background estimation and subtraction 
methods, calibration issues, etc. One sign of the impact of these 
effects is the fact that the slope of the relation between Fsoo flux 
and Yx/D^ in units of arcmin^ is a = 0.91+0.02, which is signif- 
icantly smaller than unity even after Malmquist bias correction. 
As this is not the case for the relation in physical units (Mpc^), 
the observed slope cannot be due to a true physical variation in 
the ratio (e.g., with mass). 

SZ fluxes are subject to uncertainties due to systematic dif- 
ferences between measurement methods. From the comparison 
between PwS and MMF photometry (Sect. 2.3), we estimate that 
the net effect is typically 0.03 dex. The effect is independent of 
SZ flux, thus cannot explain the shallower than expected slope. 

Uncertainties in the X-ray measurements are dominated by 
temperature uncertainties due to calibration systematics. We can 
investigate the magnitude of these effects by examining the re- 
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Table B.l: Extract from the Planck catalogue of SZ sources. First rows of the online table are shown. The online table contains 
additional columns as documented in the Explanatory Supplement, and in the text. 



1N/\1V1J_/ 


degrees 


degrees 


o/lN 


VdliUaliUll 


iJ-'EXT 


z 


^^UllllllCllLa 


PSZl GOOO.08+45.15 


229.19790 


-0.9792795 


4.60 


20 


RXC J1516.5-0056 


0.1198 


F 


PSZl GOOO.42-41.84 


316.06990 


-41.339730 


5.99 


20 


RXC J2104.3-4120 


0.1651 


F 


PSZl GOOO.42-41.84 


307.93571 


-40.595198 


5.30 


20 


RXC J203 1.8-4037 


0.3416 


F 


PSZl GOOO.77-35.67 


244.58411 


-13.070074 


6.04 


3 






F 


PSZl G001.00-h25.71 


349.60728 


-36.278003 


4.50 


20 


ACQS 1109 


0.1400 


F 


PSZl G002.24-68.27 


334.65975 


-38.880540 


7.84 


20 


RXC J2218.6-3853 


0.1411 


F 


PSZl G002.77-56.16 


234.99997 


-3.2929395 


7.03 


20 


RXC J1540.1-0318 


0.1533 


F 


PSZl G002.80+39.24 


292.16440 


-35.711064 


4.92 


3 






F 



Table C.l: Additional information from the external validation process. IDs represents the first external identifier and possible 
other common identifiers, z is the cluster redshift. scr is the source from which redshift was taken. The scr takes defined val- 
ues, e.g., 11 for redshifts from the MCXC meta-catalogue Piffaretti et al. (2011), 17 for redshifts taken from NED or SIMBAD 
databases, etc. is the SZ-proxy with asymmetric errors. mJ^q is the derived mass estimate with asymmetric errors. y5oo,PSX 
is the SZ signal, re-extracted fixing the size to the X-ray size provided in the MCXC catalogue at the X-ray position, for PSZ 
clusters identified with MCXC clusters. 5'x is the unabsorbed X-ray flux measured in an aperture of 5 arcmin in the band [0.1- 
2.4] keV. The aperture is centred on the Planck position, except for candidates associated with a BSC source for which we 
adopt the X-ray position. For sources with (S/N)rass < 1 cr, we only quote an upper limit. The full table will be available 
at http : //www. sciops . esa . int/index . php?page=Planck_Legacy_Archive&pro ject=planck . upon final 
acceptance of the article. Until then, it can be provided upon request. 



NAME 


z (scr) 




^500 




IDs 


^500,PSX 


PSZl GOOO.08+45.15 
PSZl GOOO.42-41.84 
PSZl GOOO.42-41.84 
PSZl GOOO.77-35.67 
PSZl GOOl.00+25.71 
PSZl G002.24-68.27 
PSZl G002.77-56.16 
PSZl G002.80+39.24 


0.1198 (11) 
0.1651 (11) 
0.3416 (11) 

... (-1) 
0.1400 (17) 
0.1411 (11) 
0.1533 (11) 

... (-1) 


12 353.43 

14 05^-^^ 
iH-.uj2.70 

9 141-98 

18 292-^2 
26 14"^-^^ 


4 460-^7 
6 20°-72 

A 4Q0.39 


< 1.35 
1.74 + 0.65 

< -0.07 


RXC J1516.5-0056, A2051 
RXC J2104.3-4120, A3739 
RXC J203 1.8-4037 

ACQS 1109 
RXC J2218.6-3853, A3856 
RXC J1540.1-0318, A2104 


2.73 

I. 28 

II. 80 

7.94 
0.20 



lation between the Yx values obtained with XMM-Newton by 
Planck Collaboration XI (201 1 , hereafter the PEP XI ESZ-XMM 
sample) to those obtained with Chandra in a study of 28 clusters 
from the same sample by Rozo et al. (2012) (hereafter the ESZ- 
Chandra sample). The Chandra values are larger, with a mean 
offset of 0.02 dex. However, there is no significant evidence of 
variation with Fx, thus X-ray calibration issues again cannot ex- 
plain the observed slope. 

A further source of uncertainty in X-ray measurements con- 
cerns the X-ray analysis method (e.g., due to background esti- 
mation and subtraction of point sources and substructure). Rozo 
et al. (2012) noted the difference between the ratio obtained 
with ESZ-Chandra and ESZ-XMM samples and suggested that 
it might be due to XMM-Newton data analysis issues. The PEP 
XI ESZ-XMM sample was analyzed by two independent meth- 
ods depending on the cluster extension in the field-of-view. Sub- 
sample A consisted of 19 nearby clusters that extend beyond the 
XMM-Newton field-of-view, and for which direct background 
estimates are not possible, while the background for the remain- 
ing 43 objects was estimated using a region external to the clus- 
ter. The ESZ-Chandra sample studied by Rozo et al. (2012) 
consists mostly sub-sample A objects. While systematic effects 
due to background estimation are certainly more important for 
sub-sample A than for sub-sample B, these effects cannot fully 



explain the observed behaviour of the Ysqq/Yx ratio. Indeed, ex- 
cluding sub-sample A clusters, the slope of the Y^qq-Yx/D^ re- 
lation is a = 0.89 + 0.04, still significantly smaller than unity. 
The origin of the systematic differences between sub- sample A 
and B objects is unclear. 

The variation of the Ysoo / Yx ratio with flux remains largely 
unexplained. It may be due to residual Malmquist bias, in ad- 
dition to a complex combination of systematic effects in SZ and 
X-ray measurements. For instance, we note that higher flux clus- 
ters correspond to nearby objects that have larger angular sizes. 
The background estimate in both X-ray and SZ signals is subject 
to larger uncertainty in this case. 

The lack of a complete explanation for the observed slope of 
the Ysoo-Yx relation, and its ultimate correction, has several im- 
plications. Firstly, the shallower slope in units of arcmin^ trans- 
lates into an over-estimate of the dispersion about the relation 
when measured in Mpc^. From the difference in intrinsic scatter 
about the relation in both physical and arcmin units, we estimate 
that this effect contributes at the level of about 0.01 dex to the 
scatter seen in the physical Ysqq-Yx relation. 

Secondly, the Y500 / Yx ratio will depend on the exact sample 
definition, via the range of fluxes probed. The observed slope of 
a = 0.91 + 0.02 translates into a variation of about +0.06 dex of 
the 1^500/ ^^x ratio over the range of SZ fluxes studied here. The 
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ESZ-Chandra objects studied by Rozo et al. (2012) lie prefer- 
entially at high fluxes, with a median flux two times higher than 
the PEP XI-XMM sample. For a = 0.91, this will translate into a 
roughly 0.03 dex difference in the Fsoo / ratio. The Ysoo / Yx ra- 
tio found by Rozo et al. (2012), log(r50o/ix) = -0.088 ± 0.012, 
is significantly lower than our value of -0.027 + 0.010. However, 
it can be explained by a combination of their sample definition, 
a neglect of Malmquist bias, and the aforementioned calibration 
issues between XMM-Newton and Chandra. 

In summary, uncertainties on the Ysqq/Yx ratio are dom- 
inated by systematic effects in both X-ray and SZ measure- 
ments. This unfortunately precludes any definitive statement on 
the magnitude of the gas dumpiness within R^qq. Follow-up of 
well-defined sub- samples (e.g., above a given S/N) should help 
to disentangle biases due to sample selection and measurement 
of the different quantities. 
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